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SECTION  I 


INTRODUCTION 

This  report  derives  an  integral  equation  for  the  linearized 
steady  supersonic  potential  flow  over  a  wing.  4s  in  the  case  of  a 
subsonic  flow  a  difficulty  arises  because  of  tne  singularity  of 
the  kernel  of  the  integral  equation.  The  singularity  makes  it 
necessary  to  express  the  upwash  at  the  wing  by  a  limiting  process 
in  which  one  approaches  the  wing  from  the  interior  of  the  flow 
field.  Each  integral  equation  formulation  is  based  on  the  use  of 
a  fundamental  solution.  The  idea  of  using  such  fundamental 
solutions  for  the  representation  of  the  flow  field  is  quite  old. 

It  appears  for  instance  in  so  called  "box  methods".  These  methods 
start  immediately  with  a  discretization  of  the  problem  and  in  this 
manner  reduce  the  difficulties  caused  by  the  singularity  of  the 
kernel.  In  contrast  the  present  approach  uses  analytical  means  to 
transform  the  initial  formulation  so  that  in  the  final  formulation 

one  deals  with  tractable  quantities  and  the  limiting  process  in 

* 

which  one  approaches  the  planform  no  longer  appears  explicitly.  4 
discretization  must,  of  course,  again  be  carried  out,  but  it 
appears  only  after  the  analytical  preparation  has  been  carried 
out.  The  author  believes  that  this  procedure  gives  greater 
flexibility  in  taking  the  particularities  of  a  specific  problem 
into  account.  Of  course,  one  may  then  arrive  at  a  numerical 
procedure  that  is  less  automatic. 

The  fundamental  solution  for  supersonic  flows  is  identical 
with  that  for  subsonic  flows,  except  that  certain  terms  change 
their  signs.  If  (x,y,z)  is  a  point  of  the  flow  field  and  (5,n)  a 
point  of  the  planform  then  the  fundamental  solution  is  given  by 
1/r,  where  r  *  [(x-£)2  +  (1-M2)(n_y)2  +  ( 1 -M2)z2]1 ^2.  In  the 

present  context,  x  is  greater  than  £;  the  author  has  therefore 

2  2 
consistently  written  (x-C)  instead  of  (£-x)  ) 

2 

In  subsonic  flow  the  coefficient  (1-M  )  is  positive,  in 
supersonic  flow  it  is  negative.  This  has  a  profound  effect.  In 
subsonic  flows  all  points  of  the  wing  and  of  the  wake  have  an 
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effect  on  the  potential  in  the  field,  in  supersonic  flow  only 
points  of  the  wing  and  wake  have  ail  influence  that  lie  within  the 
forecone  of  the  point  (x,y,z),  that  is  points  (£,n)  for  which 

[(x-o2  -  (M2-1)  (n-y)2-(M2-1)z2]  >  0  and  x-£  >  0 

This  is  a  significant  simplification.  A  singularity  is 
encountered  if  this  expression  is  zero.  Roughly  speaking  one 
expresses  the  potential  by  the  integral  over  5  and  n  of  such 
fundamental  solutions  multiplied  by  a  weight  function  that  depends 
upon  5  and  n.  In  a  subsonic  flow,  if  one  evaluates  the  potential 
(or  rather  its  derivatives)  at  a  point  (x,y,z)  then  r  *•  0  for 
every  point  of  the  flow  field  for  which  z  *  0.  If  the  point 
(x,yfz)  lies  on  the  planform  (in  the  £n-plane,  i.e.,  for  z  -  0), 
one  has  a  singularity  of  the  fundamental  solution  at  £  -  x,  n  «  y. 
This  makes  the  limiting  process  z  ■*  0  necessary. 

In  supersonic  flow  the  expression  r  vanishes  along  the 
entire  surface  of  the  cone  r  >  0  with  tip  at  point  x,y,z.  For 
z  A  0  the  evaluation  of  the  potential  or  its  derivatives  must  take 
into  account  the  contribution  of  the  fundamental  solutions,  with 
values  of  5,n  of  the  planform,  for  which  r  >  0.  These  values  lie 
within  the  hyperbola  in  the  Sn-plane  given  by 

(x-o2-  (m2-i ) (n-y)2  -  (M2~i)z2,  x-k  >  o. 

At  this  hyperbola  one  has  r  =  0,  the  integrand  of  the  integral 
equation  is  singular  along  this  curve,  even  for  z  A  0.  The 
singularity  becomes  more  pronounced  if  one  takes  derivatives. 

This  is  remedied  by  a  single  transformation,  but  then  the  image  of 
the  leading  edges  under  this  transformation  will  be  z-dependent. 
The  effect  is  minor  because  for  the  problem  treated  here  the 
weight  function  for  the  fundamental  solutions  is  zero  at  the 
leading  edge. 

In  Section  II  of  this  report,  the  approaches  to  the  problem 
by  the  velocity  and  by  the  acceleration  potential  are  compared 
with  each  other.  For  supersonic  flow,  at  least,  the  author  favors 
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the  velocity  potential,  because  the  pertinent  fundamental  solution 
doea  not  introduce  a  singularity  at  the  wake  of  the  points  (£,n) 
of  the  planform  for  which  it  is  defined.  Such  singularities  make 
the  upwash  field  less  smooth,  when  one  discretizes  the  integral 
equation.  This  generates  an  uncertainty,  if  one  imposes  upwash 
conditions  av  individual  control  points.  A  main  advantage  of  the 
acceleration  potential  lies  in  the  fact  that  it  anticipates  the 
contribution  of  the  wake.  But  in  most  supersonic  flows  the  wake 
has  no  influence  on  the  pressure  distribution  over  the  air  foil. 

In  Section  III  the  basic  equations  are  compiled.  This 
includes  the  introduction  of  the  Lorentz  transform  which  will 
prove  useful  in  the  treatment  of  conical  fields.  Section  IV 
describes  the  problem  in  general  terms.  In  Section  V  the  analysis 
of  the  general  case  is  prepared  by  a  two-dimensional  example.  The 
general  integral  equation  is  derived  in  Section  VI.  This 
discussion  deals  mainly  with  the  leading  concepts.  Some  detailed 
investigations  which,  although  desirable  for  mathematical 
completeness,  only  confirm  results  that  can  be  expected,  are  found 
in  Appendix  A.  Section  VII  applies  the  integral  equation  found  in 
Section  VI  to  conical  fields,  they  occur  in  the  vicinity  of  the 
tip  of  the  airfoil.  Here  the  Lorentz  transformation  is  applied  to 
suppress  an  unpleasant  numerical  singularity.  In  Section  VIII 
some  remarks  about  a  possible  discretization  are  made.  The 
treatment  of  contribution  of  points  (£,ti)  at  a  distance  from  the 
point  (x,y)  for  whicn  the  upwash  is  evaluated  is  rather  sketchy; 
more  details  are  given  for  the  contributions  of  points  (£,n)  in 
the  immediate  vicinity  of  the  point  (x,y). 

Mathematical  details  which  would  have  been  too  disruptive  in 
the  presentations  of  the  main  ideas  are  found  in  a  number  of 
appendices. 
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SECTION  II 


COMPARISON  BETWEEN  THE  VELOCITY  AND  THE  ACCELERATION  POTENTIAL 

A 3  is  well-known,  the  analytical  formulations  for  the 
linearized  flow  over  a  wing  in  terms  of  the  velocity  potential  and 
in  terms  of  the  acceleration  potential  are  based  on  the  same 
simplifications  and  are,  therefore,  equivalent.  Differences 
arise,  however,  in  the  numerical  realization.  Here  are  the  main 
points . 

1.  The  velocity  potential  arises  from  the  acceleration 
potential  by  an  integration  along  the  stream  lines  (of  the 
undisturbed  flow).  A  piecewise  linear  approximation  for  the 
velocity  potential,  therefore,  gives  a  piecewise  constant 
function  for  the  acceleration  potential  and  with  it,  at 
least  in  steady  flows,  piecewise  constant  pressures.  To 
obtain  approximations  of  the  same  quality  one,  therefore, 
needs  smoother  approximating  functions  for  the  velocity 
potential  than  for  the  acceleration  potential. 

2.  The  acceleration  potential  gives  directly  the  pressures.  If 
one  expresses  the  acceleration  potential  by  a  distribution 
of  dipoles  over  the  wing,  then  the  local  intensity  of  the 
dipoles  gives  the  pressure  difference  between  the  upper  and 
lower  sides.  In  the  wake  there  is  no  pressure  difference, 
the  unknown  dipole  distribution  is  applied  only  at  the  wing. 
In  contrast  there  is  a  jump  in  the  velocity  potential  not 
only  at  the  wing  but  also  between  the  upper  and  lower  sides 
of  the  wake  and  one  must  allow  for  doublets  over  the  surface 
of  the  wake  too.  This  is  no  insurmountable  obstacle,  but  in 
conjunction  with  the  fact  that  the  acceleration  potential' 
gives  the  desired  pressure  distribution  directly  it  makes 
the  acceleration  potential  attractive  in  subsonic  flows.  In 
a  supersonic  flow  with  a  supersonic  trailing  edge  the  wake 
has  no  influence  on  the  pressure  distribution  over  the  wing; 
but  even  with  a  subsonic  trailing  edge  only  a  small  portion 
of  the  wake  has  an  influence.  The  advantage  which  one  is 
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inclined  to  ascribe  to  the  acceleration  potential  in  a 
subsonic  flow  is  greatly  reduced  in  a  supersonic  flow. 

3.  The  integral  equations  for  both  the  acceleration  and  the 
velocity  potential  are  obtained  by  representing  the  flow 
field  by  means  of  a  dipole  distribution,  over  the  wing  for 
the  acceleration  potential  and  over  the  wing  and  wake 
surface  for  the  velocity  potential.  In  discretizing  the 
problem,  one  frequently  divides  the  surface  into  elements 
and  represents  the  dipole  density  in  the  individual  elements 
by  simple  expressions,  for  instance  by  polynomials.  At  the 
element  boundaries  there  will  be  some  discontinuity  in  the 
first  or  higher  derivatives.  The  upwash  generated  by  such 
distributions  will  reflect  these  discontinuities.  For  both 
the  acceleration  and  the  velocity  potential,  discontinuities 
will  occur  at  the  element  boundaries,  but  for  the 
acceleration  potential  one  has  in  addition  discontinuities 
along  the  boundaries  of  the  wake  pertaining  to  that  element. 
An  upwash  field  arising  in  this  manner  is,  therefore, 
traversed  by  lines  of  discontinuity  in  the  direction  of  the 
x  axis.  This  causes  an  uncertainty  if  one  uses  control 
points  to  match  the  upwash  generated  by  the  dipole 
distribution  and  tha  upwash  given  by  the  boundary 
conditions . 
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SECTION  III 


BASIC  EQUATIONS  FOR  LINEARIZED  STEADY  SUPERSONIC  FLOWS: 
LORENTZ  TRANSFORMATION 


Let  x,y,z  be  a  Cartesian  System  of  coordinates,  in  which  x 

has  the  free  stream  direction  and  the  planform  of  the  wing  lies 

in  the  x,y  plan.  Let  U  be  the  stream  velocity,  M  the  supersonic 

2  1/2 

free  stream  Mach  number,  8  -  (M  -1)  ,  L  a  characteristic  length 

and  <ji(x,y,z)  the  perturbation  potential,  which  describes  the 
deviation  of  the  velocity  field  from  a  parallel  flow.  This 
potential  is  governed  by  the  equation 


2- 

■8  <b--  +  d>--  +  d>  — 
Txx  vyy  vzz 


0 


We  . ntroduce  dimensionless  variables 


x  «  xL 

y  -  yL 

z  »  zL 

<p(x,y  ,z)  -  UL<j>(x,y,z) 


x  -  x/L 

y  -  y/L 

_  "  /♦ 

Z  -  Z/Lj 

UL<f>(x/L,y/L,  z/L) 


(1) 


(2) 


One  obtains 


“^xx  +  ^yy  ^zz  =  0  (3) 

Originally  the  author  had  introduced  at  this  point  a  Prandtl 
Glauert  coordinate  distortion  by  which  the  last  equation  is 

normalized  to  the  case  where  M  -  /2.  In  view  of  a  later 
extension  to  unsteady  flows,  where  this  distortion  loses  its 
usefulness,  this  is  not  done  here.  (It  will  be  used  in 
Appendix  A,  which  deals  with  certain  mathematical  questions,  but 
has  no  bearing  on  numerical  work.) 
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At  the  planform  the  z  component  of  the  velocity  is 
determined.  Let  the  shape  of  the  wing,  deluding  the  effect  of 
an  angle  of  attack,  be  given  by 

z  -  F(x,y) 

Then  one  has  the  boundary  conditions 

4>-(x,y ,z«0)  -  UF-(x.y) 

Hence  in  dimensionless  form 

4>z  ( x  ,y  ,z-0)  -  w(x,y)  (4) 

with 

w(x,y)  -  F-(Lx ,Ly)  (5) 

A  fundamental  solution  corresponding  to  a  source  (whose 
strength  need  not  be  defined)  at  a  point  x  -  £,  y  -  n,  z-Ois 
given  by 

<j>3(x,y,z)  -  [ (x-£) 2  -  B2(n-y)2  -  B2z2]-1/2 

provided  that  the  point  (x,y,z)  lies  in  the  aftercone  of  the 
point  (£;,  n.  0);  outside  the  afterCone  <J>S  *  0.  The  interior  of 
the  aftercone  is  given  by 

(x-O2  ~  B2(n-y)2  >  82z2  ,  xH 

The  flow  over  a  wing  with  supersonic  leading  and  trailing 
edges  can  be  described  by  a  distribution  of  such  particular 
solutions  over  the  planform 

$(s)(x,y,z)  -  ff - hU,_n)d^dn -  (6) 

A  [<x-02-B2(n-y)2-82z2]1/2 
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For  given  (x,y,z)  the  area  of  integration  is  the  portion  of  the 
planforra  which  lies  within  the  forecone  of  the  point  x,y,z.  Except 
for  the  leading  edge  the  boundary  of  the  region  is  then  given  by  the 
curve  of  intersection  of  the  forecone  with  the  £,n_plane,  i.e.,  by 

(x-£)2  -  82(n-y)2  -  82z2«  0 

This  is  a  hyperbola  in  the  5n_plane.  The  asymptotes  are  the 
straight  lines  through  the  point  5-x,  n-y  with  the  slopes 
— p  -  _+  B"1.  The  vertex  of  the  hyperbola  lies  at 

5  -  x  -$z 

n  ■  y 

For  z-0  the  hyperbola  degenerates  into  its  asymptotes. 

One  can  include  problems  with  a  subsonic  leading  edge  in  a 
procedure  based  on  Eq.  (6)  (where  the  potential  is  expressed  by  a 
distribution  of  sources)  by  introducing  a  source  distribution  on  a 
diaphragm  extending  from  the  leading  edge  to  the  Mach  wave  that 
forms  the  boundary  of  the  region  of  influence  of  the  wing.  In  the 
present  report  this  possibility  will  not  be  explored.  We  shall  use 
instead  a  representation  of  the  flow  field  by  a  superposition  of 
doublets . 

In  conjunction  with  the  treatment  of  conical  fields  (Section 
VII)  it  will  be  useful  to  carry  out  a  Lorentz  transformation* 


♦The  idea  to  use  a  Lorentz  transformation  in  the  theory  of 
linearized  supersonic  flow  is  fairly  old.  The  author  regrets  that  he 
is  unable  to  give  the  reference  in  which  it  was  first  proposed. 
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(7) 


x 

y 

X 

y 


(1-o2)-'/2x  *  6c(1-02)-,/2y 


r’cd-c2)-1'2*  .  (1  -l)-1/2j 


(1-o2)',/2x  -  ec(1-c2)-1/2y' 


■6-'c(1-o2)-1/2x’  .  (,-o2)-,/2y‘ 


z  »  z 

4>(x,y,z)-4>(x,y,z) 


(8) 


The  coordinates  £  and  n  are  transformed  in  the  same  manner. 
Then 


-B2* 


xx 


yy 


zz 


2- 

“6  4)  +  d>  +  A 

rxx  ^yy  Tzz 


The  hyperbola 


(x-£)2  -g2(n-y)2  -B2z2 


transforms  into 

(x-£)2  -  &2(n-y)2  -  B2z2  -  0  (9) 

and  one  obtains  for  the  Jacobian  of  the  transformation 


8<x,y) 

After  the  Lorentz  transformation  one  has 
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(10) 


-  ff - - 

Cx-5)2  -  B2( n-y)2  ~  B2z2]1/2 
with  hU.n)  -  hUU,n)  ,nU,n)) 

The  functions  5(5, n)  and  n(5,n)  are  given  by  expression  analogous 
to  Equations  (8).  The  boundary  of  the  region  of  integration 
formed  by  the  hyperbola,  appears  in  the  same  form  as  before, 
i.e.,  one  obtains  Eq.  (9).  The  boundary  formed  by  the  leading 
edge  must  be  expressed  in  terms  of  the  coordinates  £  and  n. 

If,  for  instance,  the  leading  edge  is  given  by 

y  -  c^  x  -  0,  with  c1 >B  1 

then  one  obtains  in  terms  of  x  and  y 

y(1 +Bcc1 )  -  x(c1  +  B_1c) 


Setting 


c 


1 


one  obtains  as  equation  of  the  leading  edge 

x  *  0 

This  is  a  significant  simplification. 
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SECTION  IV 


GENERAL  REMARKS 

For  a  wing  of  zero  thickness  and  for  boundary  conditions 
satisfied  at  the  planform,  the  prescribed  z  component  of  the 
velocity  is  symmetric  with  respect  to  the  plane  z  -  0  (this  is 
the  plane  of  the  planform).  The  potential  is  then  antisymmetric. 
Outside  the  wing  and  the  wake  the  potential  is,  therefore,  zero; 
within  the  wing  and  the  wake  a  jump  of  the  potential  from  some 
value  to  the  value  with  the  opposite  sign  must  be  admitted.  The 
pressure  distribution  over  the  wing  is  influenced  only  by  those 
points  of  the  wake  for  which  part  of  the  aftercone  lies  within 
the  wing  surface;  because  the  effect  of  the  upwash  at  any  point 
of  the  flow  field  is  felt  only  in  its  aftercone.  For  a 
supersonic  trailing  edge  the  wake  has  no  effect  on  the  wing.  A 
potential  which  has  the  desired  antisymmetry  is  obtained  by  a 
distribution  of  doublets  (oriented  in  the  z-direction)  over  the 
wing  and,  if  necessary  over  part  of  the  wake.  A  doublet 
potential  is  obtained  by  a  differentiation  of  the  source 
potential  with  respect  to  z.  In  formulating  the  boundary 
conditions  at  the  wing  one  must  express  the  z-component  of  the 
velocity  field.  One  thus  encounters  a  first  differentiation  with 
respect  to  z  when  one  derives  a  doublet  potential  from  the  source 
potential  and  a  second  differentiation  when  one  formulates  the 
boundary  conditions.  The  singularities  which  arise  by  these 
differentiations  are  the  main  concern  of  the  following 
discussions.  In  the  subsonic  case  these  singularities  occur  only 
for  5  «  x  and  n  *  y.  In  contrast  they  occur  in  the  supersonic 
case  also  along  the  hyperbola  which  form  the  boundary  of  the 
region  of  integration. 

In  the  subsonic  case  one  has  the  corresponding  potential  for 
a  source 


with 


1 1 


2  p  2  2  2  1  /2  2  2 

r  -[(x-5)  +  B  (n-y)  +  B  z  ]  with  B  -  l-M* 

Here  r  represents,  in  a  suitable  metric,  the  distance  between  the 
point  (x,y,z)  and  the  point  (C,n,0).  This  can  be  carried  over  to 
the  supersonic  case  if  one  introduces  as  distance  definition 

[(x-O2  -  B2(y-n)2  -  z2]1/2 

The  fundamental  solutions  are  singular  at  all  points  where  this 
"hyperbolic"  distance  is  zero  and  this  includes  the  points  of  the 
bounding  hyperbola. 

Incidentally,  in  a  field  with  hyperbolic  distance  definition 
the  Lorentz  transformation  is  the  counterpart  to  a  rotation  in  a 
field  with  the  elliptic  distance  definition. 
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SECTION  V 


THE  TWO-DIMENSIONAL  CASE 


To  familiarize  ourselves  with  the  mathematical  technique  we 
first  treat  the  two-dimensional  case.  If  one  introduces  the 
Lcrentz  transformation  this  includes  the  treatment  of  the 
infinitely  long  swept  wing  with  a  supersonic  leading  edge.  For 
the  two-dimensional  case,  the  leading  edge  is  always  supersonic. 
The  function  h  does  not  depend  upon  y  and  the  evaluation  of  the 
potential  can  be  carried  out  for  y  -  0.  The  upwash  is  needed  for 
z  -  0.  One  has 


♦  <  x ,z) 


h(£)d£dn 

1 


(id 


As  wa3  stated  above  the  area  of  integration  is  bounded  by  the 
leading  edge  (here  taken  as  the  line  £  -  0)  and  the  intersection 
of  the  forecone  of  the  point  (x,y,z)  with  the  plane.  Because 
we  carry  out  the  integration  for  y  -  0,  the  bounding  hyperbola  is 
given  by 

( X- C) 2  -  82n 2  -  B2*2  -  0 

In  the  limit  z  -  0  the  boundary  is  given  by  its  asymptotes 

x  -  E;  +  Bn  -  0 
The  vertex  of  the  hyperbola  lies  at 

x  -  £  ■  Bz ,  n  °  0 

With  the  limits  specified,  Eq.  (11)  appears  in  the  form 


$(x 


x-8z  upper 

,0,z)«*jh(£)(j  - 


(x-O 


e2z2  -  e2ri2]1/2 


lower 


where 


r/  2  .2  2 -,1  /2  /D 

Slower  =  "  [(x'^  “  6  z  ]  /e 
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r ,  2  .2  2-.1  /2  ,D 
"upper  ’  *  C<x-5)  ■  B  z  3  /B 


Setting 


one  obtains 


-  -1r/  rN2  .2  2-.1/2 
n  -  riB  [(x-O  -B  z  ] 


:u2-v2  -  s2*2  -  bV]  -  [<x2-o2  - 

and  for  the  inner  integral 


+  1 

»-l[  p"1-^ 


x-Bz 


Hence  $(x 


|Z)  -  B_1ir| 


h(£)d£ 


<j>x  -  B~  ir  h(x-Bz) 

<f>  »  -  irh(x-Bz) 

Z 

This  satisfies  the  differential  equation  for  one- dimensional  flow 


~  8  (b  +  A  =  0 
yxx  yzz 
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SECTION  VI 


THE  UPWASH  IN  THE  GENERAL  CASE 

The  point  of  departure  is  again  Eq.  (6).  Let  the  leading 
edge  be  given  by 


n  -  g(5)  (12) 

We  set,  in  analogy  to  the  procedure  of  the  two-dimensional 
problems 


n  **  y  +  ne_1q(x,£,z)  (13) 

n  -  B(n-y)/q 


where 

q(x-g»s)  -  [(x-5)2-e23]1/2  (U) 

To  emphasize  that  in  the  expression  q  the  coordinate  z  occurs 

2  2 
only  in  the  form  of  z  ,  we  have  introduced  z  »s 

One  has  q(x-£,0)  =  (x-O 

H -(-«'«  (,5) 

If  •  -82*/q 

The  portion  of  the  boundary  of  the  region  of  integration 
formed  by  the  hyperbola  is  transformed  into 


T|  »  +  1 


15 


The  portion  of  the  boundary  formed  by  the  leading  edge  appears  as 


n  -  g(xfy,CiS) 


(16) 


where 


g(x,y,£,s) 


B(g(P-y) 

qCx-T.3  J 


(17) 


The  expression  for  generated  by  a  source  distribution  is  then 
given  by 

x-8z  nupper 

(f  dnut  ns) 

l  (1  -h2)1/2 

nlower 


<j>(s)(x  ,y  ,z)  -  8~ 1 1 

5o 


Here  £q  is  the  smallest  value  of  £  in  the  region  of  integration. 
If  the  wing  has  a  tip,  £Q  is  independent  of  z.  In  cases  where 
the  foremost  point  of  the  region  of  integration  is  the 
intersection  of  the  hyperbola  with  the  leading  edge,  £0  will 
depend  upon  z.  Furthermore, 


nlow*r  “  +1  alon8  hyperbola  (19) 

upper 


n, *  g(x,y,£,z)  if  the  limit  of  n  lies  at  the 

jl  ow  cr 

upper 

leading  edge. 


In  Eq.  (18)  the  limit  z+0  can  be  formed  immediately 


x  nupper 

( |  dn)d£  (20) 

^lower 


<|>^(x,y,z-0) 


1 


5o 
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This  formulation  is  practical  for  supersonic  leading  and  trailing 
edges.  At  the  planform  the  upwash  is  given  by  Eq.  (4). 
Accordingly,  one  must  form  in  Eq.  (18)  the  derivative 
4 ,y tZ“0) .  Before  we  do  this  let  us  determine  which  is 

needed  to  express  the  pressures  in  cases  where  one  represents  the 
flow  field  by  a  superposition  of  sources.  If  one  is  interested 
in  only  the  pressures  of  the  planform,  then  one  can  set  z-0 
before  one  carries  out  the  differentiation,  i.e.,  one  can 
differentiate  Eq.  (20).  The  derivative  with  respect  to  the  upper 
limit  of  the  outer  integral  gives 


B_1h(x  ,y ) 


dn 


(1 


■'n2>,/2 


B  1ivh(x,y) 


We  denote  by  h^  and  h^2^  the  partial  derivatives  of  h  with 
respect  to  its  first  and  second  arguments.  One  then  obtain 


x  ^upper 

<f»®(x,y,z-0)  -  B_1h(x,y)  ♦  B"2J  ( J ^ (21) 

n 

lower 


The  limits  for  n  are  determined  for  z-0.  Hence  from  Eqs.  (19) 
and  (13)  and  (14) 


nlower  ’ 
upper 


or 


T|  lower 
upper 


pg(U-y 
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( s ) 

Next  the  derivatives  of  $  with  respect  to  z  are 
evaluated.  As  we  mentioned  above,  the  first  derivative  is  needed 
in  satisfying  the  upwash  conditions,  if  one  expresses  the 

(  3  ) 

potential  by  a  source  distribution.  Besides,  the  expression  4>z 
can  be  interpreted  as  the  potential  due  to  a  doublet 
distribution. 


d  s 

<f>  (x,y,z)  -  <j>z(x ,y ,z ) 

To  formulate  the  upwash  conditions  for  the  potential  <j>d  one  must 
form 


<J>zU,y,z)  -  <|»|z(x,y,z) 

Ultimately,  all  quantities  will  be  evaluated  for  z-0. 
Configurations  of  the  region  of  integration  are  shown  in  Figures 
1  and  2.  The  leading  edge  may  contain  subsonic  or  supersonic 
parts.  For  supersonic  or  subsonic  parts  of  the  leading  edge  the 
slope  dn/d£  *  dg/d£  is  respectively  greater  or  smaller  than  B  1 • 
A  leading  edge  which  is  entirely  subsonic  consists  of  two  parts, 
one  with  positive  the  other  with  negative  slope.  We  divide  the 
region  of  integration  into  two  parts  I  and  II  by  a  straight  line 
lying  in  the  region  where  the  upper  and  lower  limits  of  n 
are  +J  .  The  contributions  to  the  potential  from  these  regions 
are  denoted  by  and  i}/3'11^  or  <|>(d,I)  and  <j/d,I3:',  (if  they 

are  generated  respectively  by  sources  or  doublets),  Figure  3. 

Consider 
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Then  by  differentiation  with  respect  to  z 


4>IS,I)(x,y,z) 


i  -  - 1 

fh(E,y+nB  Q)dn 

j - — 

-1  (1  n  *  E«x-B z 


x-Bz  +1  m 

■2 [  r  hv  '(E.y+nB  g)  n(3g/3z)dnjdE 

J  ( 1  - n2 ) 1  /2 

*1 


Because 


q(x-E.z) I  -  0 

E  -  x-Bz 


the  first  term  gives  immediately 


-h(x-Bz 


+  1 

,y)  f - r 

(1-n 


2-f/2  *  "  irh(x-Bz.y) 


Substituting  from  Eq.  (15)  one  obtains  for  the  second  term 


x-Bz  +1  (0s  , 

"I  <1 

51  "I 


The  factor  q  tends  to  zero  as  E  approaches  its  upper  limit 

and  for  z«0  (s-0)  one  obtains  a  non-convergent  integral  of  the 

form  l^TF'  One  observes,  however,  that 
X"  4 


(  h(2)U,y)ridn 


.1 


because  of  the  antisymmetry  of  the  integrand  with  respect  to  n. 
The  inner  integral  can,  therefore,  be  replaced  by 
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[  Lh^ 2^  (5 ,  (y+nB~  'q )  )-hv  (5 , 

(1-W72 


)  3nd~i 


h(2)(5,y+Q-  h(2)(5,y) 


I^Uty.C) 


The  function  exists  even  for  5-0.  There 

^U.y.O)  -  h(2,2)(5,y) 

(2  2) 

(where  h  ’  denotes  the  second  partial  derivative  with  respect 
to  the  second  argument  of  h.) 


Then 


h(2)(5,y*n8~1q)  -  h(2)U,y)  -  nB^q^  (5,y  ,nB_1q) 


(s  I) 

and  the  expression  <(r  ’  assumes  the  form 


(3,1) 


-irh(x-Bz  ,y  )-zB_1  j  (J 

Ci  - 


i"fZ  rf  ♦i^»y+qB",|q)n2dn) 


(1-n2)1/2 


The  inner  integral  is  bounded.  In  the  limit  z-0  the  second  term 
vanishes  and  one  obtains 

♦  (x,y  ,z-0)  -  <f>i3,I*(x,y  ,z=0)  -  -irh(x,y)  (23) 


The  contribution  of  the  region  I  to  the  second  derivative  of 

(  3  ) 

<f>  with  respect  to  z  is 
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4>id,I)(x,y,z)  -  irf5h(1)(x-f}z,y) 


x-Bz  + 

I1 1  < 


4^(5. y.nS  1  q ) )  ti2cIti) 

- 77"-?rf/5 - 

( l-n  ) 


)d£ 


■zB 


x-Bz 

-2  f 

J 

*1 


+1 


<!>j3)U,y,n8' 


^Jr^Uq/Sz)  _ 


77-W75' 


dn)d5 


( 1) 

Here  denotes  the  partial  derivative  of  with  respect  to 

its  third  argument.  The  second  term  on  the  right  has  already 
been  evaluated,  it  gives  a  non-vanishing  contribution.  The  third 
term  vanishes  in  the  limit  z-0,  because  of  the  factor  z.  The 
fourth  term  has  actually  a  factor  z2,  because  *  ~B2z/q.  Thus, 
it  will  vanish  in  the  limit  z+0  provided  that  the  remaining 
expression  remains  finite.  One  has  after  substitution  of  3q/3z 


x-Bz 

:  f 
J 


^|3^(5,y,q6_1q))n3dri 


-1 


( 1  -  n  ) 


-T7f72' 


Again  one  is  concerned  about  the  factor  q  in  the  denominator 
which  vanishes  at  the  upper  limit  C-x-Bz.  Here  the  procedure 
shown  above  is  applied  again.  One  observes  that 


+  j  U,y  ,0)n3dn 

!1  (T??72 


0 
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because  of  the  antisymmetry  of  the  integrand  with  respect  to 
n-0 . 

One  introduces 


^3)  U,y  U.n.O) 

^>2U, y,S)  -  - 1 - 


( 3) 

If  is  sufficiently  smooth,  as  it  is  except  for  possible 

singular  lines,  the  function  exists  for  5=0. 

Therefore, 


+  1  ^n^dn 

?!  oW75 

-1 


and  this  expression  exists. 
One  then  obtains 


+  1  _ 1 

_1  q)ndn 

s  J 


-1 


(1-T1  ) 


-2.1/2 


♦id,I)(x,y,z-0) 


*  *}  n^U.y.nB^x-OMn 

ir0h(1  ^(x.y)-^”1  J  (J  ^2U2  )d5 

£.]  -1  ( 1 -n  ) 


Substituting  ip ^ ,  one  obtains  for  the  second  term 


A 

i  ' 


-  (?)  -  - 1 

nhv  U,y+n8  (x-O)dn 

(r?p^ 


)dC 


(24) 


As  5  approaches  the  upper  limit  x  the  integrand  seems  to  tend  to 
infinity  but  by  the  above  argument  one  obtains  a  finite 
expression  for  the  inner  integral  provided  that  the  integrations 
are  carried  out  in  the  sequence  shown  in  Eq.  (24).  To  allow  for 
a  change  in  the  sequence  of  integrations  we  write 
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wvk  lorwtvn 


Now  it 


x-e  +1 


lira  f  (f  ijh^  (g,  y+n8  1  (x-&) )dnjd^ 


-„-l 


£•*■0  j  '‘j 


5l  ^ 


(1-n2)1/r(x-c) 


is  permissible  to  change  the  sequence  of  integrations 


+  1  x-e 

e.O  1  — ^TTST  (I  h<a«.(^;r.?.te-SUflS)d-w 

e  0  -1  (1-n  )’  /d  e1  5 


lira 


(25) 


According  to  this  formulation  one  first  excludes  the  neighborhood 
of  the  point  (6,n)  *  (x,y)  by  a  straight  line  5»x-e  (e>0),  and 
after  the  integrations  forms  the  limit  e+0.  Then 

♦id,I)(x,yfz-0)  -  ir$h(1)(x,y)  +  I  (26) 

where  I  is  evaluated  either  in  the  form  of  Eq.  (24)  or  of 

Eq.  (25).  The  regions  of  integrations  for  z-0  and  z*0  are  shown 

in  Figure  4. 

In  the  region  II,  x-£>0  therefore,  q-[ (x-£)  ^-B^z2]1  /2>0 , 

even  for  z-0.  A  denominator  q  is  no  longer  detrimental.  But 

here  part  of  the  boundary  of  the  region  depends,  after  the 

2 

transformation  to  n,  upon  s-z  .  (The  transformation  to  n  is 
desirable  because  of  those  portions  of  the  boundary  formed  by  the 
hyperbola.)  We  interchange  the  sequence  of  integrations  in 
Eq.  (18): 


♦(s,II)(x,y,z)  -  8_ 1  j  - ~\~XJ2  (}  h(C,y+ng_1q)d5)dn  (27) 

1  f(x,y,n,s) 

Here  f(x,y,n,s)  denotes  the  inverse  of  the  function  g(x,y,£,s), 
Eq.  (17),  at  fixed  x,y  and  s.  Then 
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*i3,II)(x,y,z)  -  s'1 


77^T7J(-Hrf  h<5.y*ns"1q)  j  )dn 


5«f(x,y,n,s)  (28) 


+  1 


f  n  rl1h<2)(5,y*Ss"1q)()E1 
'  f(x,y,n,s) 


dn 


In  the  first  expression  h  is  evaluated,  for  the  particular  value 
of  n  under  consideration,  at  the  lower  limit  of  the  region  of 
integration  over  5.  This  is  a  point  of  the  leading  edge  and 
there  h«0.  For  z»0  both  terms  on  the  right  of  the  last  equation 
vanish 


$2S,II)<x,y,z-0)  -  0 


We  combine  this  result  with  Eq.  (23)  and  obtain 

<t>|(x,y ,z«0  +  )  -  <j>d(x ,y ,z=0) *  -irh(x,y)  (29) 

( s ) 

We  now  turn  to  the  second  derivatives  of  <J>  with  respect  to  z. 
At  the  leading  edge  the  potential  and  with  it  h  is  zero  and 
continuous.  The  first  term  on  the  right  of  Eq .  (28)  is, 
therefore,  zero,  even  for  z^O.  Differentiating  the  second  term 
with  respect  to  z,  one  obtains  an  expression  which  will  not 
vanish  for  z-0,  namely 


+1 


f 

-  J 

-1 


h(2}  U,y+nP  Vx-C.O) 


q(x- £  >0) 
f (x,y,n,0) 


dfOdn 


(30) 


Further  terras  arise  by  differentiations  of  s  with  respect  to  z; 

they  introduce  other  factors  z.  One  will  surmise,  that  these 

terms  vanish  for  z=0.  This,  however,  is  not  self-evident.  At  a 

1  /2 

subsonic  leading  edge  h  behaves  as  u  ,  where  u  is  the  distance 
from  the  leading  edge.  The  differentiation  with  respect  to  z 
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generates  some  terms  which  are  infinite  as  5  approaches  the 

( 2) 

leading  edge.  (One  of  these  is  h  itself.)  Actually,  they 
will  cancel.  A  detailed  discussion  is  carried  out  in  Appendix  A. 
Eq  (30)  is  indeed  the  only  contribution  to  ^  (x  ,y  ,z=0) . 

Cm 

Anticipating  these  results,  we  compile  here  the  formulae  needed 
for  a  numerical  approach. 


4>d(x,y,0  +  )  *  -irh(x,y) 


♦f(x.y.O)  -  *Bh(,)(x,y)-!f  dfrli  (31) 

2  ]}  (  1-T1  )  /2(X“0 


iJ 

A' 


where  A'  is  the  area  of  the  planform  within  the  fore-cone  of  the 
point  (x,y,z«0),  but  with  the  point  (fj,n)-(x,y)  excluded  by  a 
line  £>x-e,  e>0.  Subsequently,  one  forms  the  limit  e+0.  The 
same  expression  is  written  in  the  original  coordinates.  Equating 
the  resulting  expression  with  the  given  downwash  w(x,y),  one 
obtains  the  integral  equation  from  which  the  doublet  density 
h(x,y)  can  be  determined.  The  limiting  process  z-»0  no  longer 
appears  in  this  equation. 


w(x,y)-4>^(x,y  ,0)-irgh(1  *  (x,y)- 


(2) 

(n-y)hv  'U , n ) d cd n 


2  r  f 

’  {i[x-U2-B2(n-y)2]1/2(x-02 


(32) 


If  the  integration  with  respect  to  n  is  carried  out  first,  then 
the  exclusion  of  the  point  (£,ti)  =  (x,y)  is  not  necessary.  But 
if  the  integrations  are  carried  out  in  a  different  sequence  then 
in  some  intermediate  steps  one  obtains  contributions  which  are 
large  if  £  approaches  x  and  the  precautionary  step  of  introducing 
A',  which  initially  excludes  the  critical  vicinity  of  the  point 
(x ,y ) ,  Is  necessary. 

Eq.  (31)  expresses  the  upwasn  in  terras  of  h^1^  -  hx  at  the 
point  (x,y)  and  of  an  integral  which  is  due  to  the  deviation  from 
a  two-dimensional  flow.  (One  notices  that  the  integrand  contains 
a  factor  hv  =  h^.)  Because  of  the  denominator 
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2  p  2  1/2 

[(x-0“  "  8  (r,-y)  ]  '  ,  points  in  the  vicinity  of  the  Mach  waves 
through  the  point  (x,y)  enter  with  greater  weight.  A  good 
approximation  for  the  vicinity  of  these  lines  is,  therefore, 
desirable . 

Some  thoughts  about  the  application  of  Eq .  (32)  in  a 
numerical  procedure  are  found  in  Section  VIII. 

All  transformation  carried  out  in  this  Section  could  have 
been  performed  in  Eq.  (10).  Therefore,  one  obtains 
alternatively 


w(x ,y )  *  <i>!?(x,y,0)  =  irgh~(x,y)-  B2[ 


(n-y)0(h(£,n)/dn)dsdn 
- - '*■ — » - j-r  j* - *r 

[(*-?)  -8^(11-?)  3  (x-C) 


h(x,y)  has  been  defined  in  conjunction  with  Eq.  (10). 


w(x,y)  -  w(x(x  ,y )  ,y(x  ,y) ) . 


Moreover 


SECTION  VII 


CONICAL  FIELDS 

Frequently,  the  planform  has  a  tip  at  which  two  straight 
subsonic  leading  edges  meet.  In  the  vicinity  of  this  tip  one 
obtains  a  conical  field.  The  original  concept  of  conical  fields 
is  applicable  if  the  angle  of  attack  is  constant.  The  wing  is 
then  regarded  as  a  degenerate  cone  generated  by  straight  lines 
through  the  tip  lying  in  a  plane.  In  this  case  the  velocity 
vector  is  constant  along  each  straight  line  through  the  tip. 
Basically,  one  deals  with  a  similarity  solution,  which  reduces 
the  number  of  dimensions  of  the  problem  from  three  to  two. 

In  linearised  flow  the  concept  of  a  conical  field  is 
applicable  also  if  the  upwash  behaves  like  a  power  of  the 
distance  from  the  tip.  Separately  for  each  such  power,  one  makes 
a  similarity  hypothesis.  If  the  velocities  are  constant  along 
straight  lines  through  the  tip  one  can  apply  a  further 
transformation  to  the  resulting  partial  differential  equations  in 
two  dimensions  (for  the  individual  velocity  components)  which 
yields  the  Laplace  equation.  In  this  case  the  theory  of  analytic 
functions  can  be  applied.  Unfortunately,  one  is  led,  at  least  in 
the  present  problem,  to  the  evaluation  of  elliptic  integrals. 
Moreover,  if  the  velocity  is  not  constant  but  follows  a  general 
power  law,  then  the  theory  of  analytic  functions  is  not 
applicable.  We  show  here,  how  such  conical  fields  can  be  treated 
(numerically)  by  means  of  the  present  integral  equation 
formulation. 

The  procedure  is  first  described  in  general  terms.  Let  the 
leading  edges  be  given  by 

y/x  -  a1  and  y/x  =  -a2  ;  a1  >  0,  a2  >  0 

The  planform  is  nearly  always  symmetric  with  respect  to  the  x 

-1  -1 

axis.  Then  a^  ■  a2.  For  subsonic  edges  a1  <  3  ,  a2  <  B 
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The  present  approach  assumes  that  the  upwash  is  given  by  a 
linear  combination  of  expressions 

w(x,y)  -  xm  w(m)(y/x)  (33) 

In  practice  one  usually  has 

w(x,y)  «  xmy^ 

Then  m  -  m+j 

ft(y/x)  -  (y/x)j  (34) 

Ttye  computations  are  carried  out  separately  for  each  choice  of  m 
and  j . 

For  a  constant  angle  of  attack  m  =  0  and  is  constant. 

In  most  practical  applications  m  =  0,  1  or  2  will  be  sufficient. 

The  unknown  function  h  is  now  written  in  the  form 

h(x,y)  =  Z  b  hmn(y/x)  (35) 

n 

where 

h,n  -  xm*V(y/x) 

The  functions  hn(y/x)  are  known.  Probably,  one  will  use  the  same 
functions  hn  for  all  values  of  m.  The  coefficients  bn  are  to  be 
determined.  Substituting  Eq .  (35)  into  the  integral  equation, 

Eq.  (32) ,  one  obtains 

xm»(m)(y/x)  =  lb  Qran(y/x)  (36) 

n 

where  Qmn(y/x)  is  the  upwash  generated  by  a  function 
hmn=xmhn(y/x) .  This  substitution  gives 

Qmn(y,x)  =  Qran1  ♦  Qmn2 
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with 


Qmn1  -  Bir(3hmn(x,y)/3x) 


,mn2 


-B 


2 f  f  (n-y)  (3hmnU,  n)/3n)d 5d _ 

JJ[(x-j-)2-e2(n-y  )23l75Tx-e) 


2 


(37) 


In  discretizing  the  problem  one  uses  only  a  finite  number  of 
functions  hn(y/x).  To  determine  the  coefficients  bn  a 
collocation  method  is  probably  sufficient.  We  set  for  this 
purpose  y/x  «  ck  and  admits  as  many  values  of  ck  as  one  uses 
values  of  n.  Then  one  obtains  the  following  system  of  equations 


xmfl(ra)(ck) 


N 

E 

1 


bnQ  <V 


(38) 


The  main  task  is  the  determination  of  the  matrix  elements 

Qran(ck). 

Regarding  the  cnoice  of  the  functions  hmn  we  make  the 
following  observation.  If  the  angle  at  the  tip  is  small,  and 
a1*a2“a,  then  one  obtains  from  slender  body  theory,  that  the 
potential,  and  with  it,  h  is  given  by 

2  2  2  1  /2 

const  (ax  -  y  )  ;  y  <  ax, 


This  suggests  that  in  a  more  general  situation  h  can  be 
represented  by 


h  =  Z  bn  hmn(y/x)  *  E  bnxm+1hn(y/x)  (39) 

with 

hn(y/x)  -  [  (a^  (y/x)  (a2+(y/x) )  31 /2  (y/x)n  (40) 

Since  for  a  slender  tip  the  first  term  is  already  sufficient 
it  can  be  expected  that  the  number  of  terms  necessary  to  obtain  a 
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reasonable  accuracy  is  fairly  small.  Moreover,  if  a1  “  a2  and 
the  upwash  is  either  symmetric  or  antisymmetric  with  respect  to 
the  x  axis,  then  only  the  values  of  n  which  respectively  are  even 
or  odd  will  occur. 

The  leading  edges  are  given  by 

y/x  -  a^  ;  0  <  a1  <  B_1 


and 


y/x  -  -a2  ;  0  <  a2  <  B 

The  further  development  is  suggested  by  the  following 
observation.  The  expression  Qmn1  defined  in  Eq.  (37),  for  m  *  0 
and  for 


(a2*2 


-  y2)1/2 


contains  an  essential  term 

3h/3x  =■  a2(a2  -  (y/x)2)-1  /2 

This  expression  and,  therefore,  Qmn1  tends  to  infinity  as  y/x 
tends  to  ....  But  one  expects  from  slender  body  theory  that  the 
upwash  caused  by  this  expression  is  finite  over  the  entire  width 
of  the  airfoil.  The  expression  Qran2  will,  therefore,  contain 
some  term  which  gives  an  equally  large  contribution  but  with  the 
opposite  sign,  so  that  this  infinity  is  cancelled.  This 
complication  can  be  avoided  if  one  first  carries  out  a  Lorentz 
transformation  (separately  for  each  of  the  chosen  values  y/x)  so 
that  the  chosen  value  of  y/x  is  transformed  into  y/x  «  0. 

The  Lorentz  transformation  is  first  carried  out  without  this 
specialization  in  order  to  provide  a  clear  picture  of  the 
individual  steps.  We  repeat  the  transformation  formulae 
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(Mia) 


x  -  (1-c2r1/2  (  x  -  ecy) 
y  -  (1-c2)~1/2  (-B_1cx  +  y) 
x  -  (1-c2)"1/2  (  x  +  Bey)  (41b) 

y  -  ( 1-c2)-1 /2  (B_1cx  +  y) 


The  variables  §  and  n  are  transformed  in  the  same  manner. 
By  the  Lorentz  transformation  the  left  hand  side  of  Eq.  (38) 
assumes  the  form 

1 -Be (y/x) 

The  constant  c  will  be  chosen  in  such  a  manner,  that,  y  *  0  for 
y/x  -  ck<  From  the  second  of  Eq.  (41b)  one  then  obtains 

c  -  -Bck  (42) 

With  this  choice  one  obtains  for  the  left  hand  side  of  Eq.  (38) 

(1-c2)_m/2  xm«m(ck)  (43) 

To  transform  the  right  hand  side  of  Eq.  (38)  one  must  first 
carry  out  the  transformation  to  £  and  n  in 

hmn(x  ,y )  -  xra+1[(ar(y/x))(a2+(y/x))]1/2(y/x)n  ;  (44) 

hmn(5,n)  is  obtained  by  substituting  x  *>  x(x,y)  and  y  -  y(x,y). 

Let  Bn/£  ■  P  (45) 


Then,  from  Eq.  (41a) 


n/£ 


P-c 

Bfl-cp) 


(46) 
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Moreover 


S  -  (1-c2)~1/2£(1~cp) 

The  leading  edges  are  originally  given  by 
n/5  -  a1  and  n/5  -  -a2 


One  finds  from  Eq. 

n/5 

The  leading  edges, 


(41) 

g  ^(y/x) 

1 +gcy/x  ' 

therefore,  transform  into 


p-a^  and  p  -  -a2 


with 


(47) 


(48) 


(49) 


Ba1 +c  ~  ga2-c 

a1  “  T+cga^  ’  a2  “  f+cga"j 

In  the  expression  hn  Eq.  (40),  a  factor 

[(a1-n/£)(a2+n/£)]1/2 


(50) 


occurs.  One  find3  with  Eq.  (45) 


1 +  ga.|  C 

(arn/5)  *  sTf-cp) 

(a^p) 

(51 ) 

1-ga0c 

(a2  n/^  ’  g(  1-cp) 

(a2+p) 

(52) 

mn 

Now  the  transformation  of  w  *  ,  Eq .  (44),  can  be  carried  out 
hranU,n)  =  (1-c2)"m/2  a3  B“n“\ra+1  Fmn(p)  (53) 
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with 


a3  -  (1-c2)“1/2[(1+ga1c)(1-Ba2c)]1/2  (54) 

Fmn(p)  -  C(arp)(a2+p)31/2(1-cp)m”n(p-c)n  (55) 

where  p  is  expressed  in  terras  of  n/£  by  Eq,  (45).  He  note 

dFmn/dp  -  [(a1-p)(a2+p)]~1 /2  fmn(p)  (56) 

where 

fmn(p)  -  (-p^^)(1-cp)m"n(p-c)n  (57) 

+  (a1  -p)  (a2~p)  [  (-c)  (m-n)  ( 1~cp)ra  n  ^(p-c)11 
+  n( 1-cp)m-n(p-c)n"1  ] 

The  second  term  in  the  last  bracket  vanishes  for  n-0,  therefore 
no  difficulty  arises,  even  for  p-c. 

One  has  in  particular 

fran(0)-(ara2)(-c)n  +  a,  a2[ (m-n)  (-c)n+1  +n(-c)n-1  ]  (58) 

Moreover 

(dfmn/dp)|  -  -(-c)n  +  (3/2)(a1-a2)[(ra-n)(-c)n+1+n(-c)n“1] 

p-0 

+  a.j a2[  (ra-n- 1 )  (m-n)  (-c)n+2+2(m-n)n(-c )n 
+  n  (n-1)(-c)n"2]  (59) 
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One  thus  obtains  for  general  c,  by  substituting  Eq.  (53) 

Qmn(x ,y )  -  (1-c2)“m/2  a3  {8"nTt3(xm+1Fmn(ey/x))/3x) 
n-n+l f  f(n-y)(3(era+1Fran(8n/5))/3n 

'  1 1  coTl)  1  x 

To  make  the  dependence  upon  5,  n  and  x,  y  obvious  we  have  here 
substituted  the  expression  for  p.  With  the  special  choice  made 
above,  namely  c  -  -0ck,  one  obtains  y  -  0  for  y/x  -  ck.  Then  the 
expression  simplifies  to 

Qran(x,0)  -  (1-c2)“m/2xmQran  (60a) 

with 

qmn  .  5mn1  ♦  Qmn2  (60b) 


where 


J“n'  -  a3$"nx"m3(im+1Fmn(0))/3x 


-a3B 


-n+1 ~-m  f  f  n3(^m''1Fmn(Bn/€))/3n)dCdn 
^  [  x-O2  -"T^n2]1  /2(x-i)2 


(61) 

(62) 


Substituting  Eqs.  (43)  and  (60)  into  Eq.  (38)  one  obtains 

wm(ck)  -  l  bnQmn(Gk)  {63) 

In  Qmn^ ,  Eq.  (61),  Fmn(0)  from  Eq.  (55)  is  substituted  and  the 
differentiation  with  respect  to  £  is  carried  out. 


Then 

Qmn1  -  a30_n(m+1)(a1a2)1/2(-c)n  (64) 
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A 

In  the  expression  Qm  ,  Eq.  (62),  p  defined  in  Eq.  (M5),  is 
introduced  as  variable  of  integration  instead  of  n.  Then 

dn  -  3  ^dp  (65) 


In  addition,  we  set 


5/x  -  £.  (66) 

The  region  of  integration  is  a  quadrangle.  Two  sides  are  formed 
by  the  leading  edges.  The  lines  p  -  const  are  straight  lines 
through  the  tip  of  the  airfoil.  The  limits  for  p  are  then 
according  to  Eqs.  (M9)  p  -  -a2  and  p  -  a1  . 

The  other  two  sides  are  given  by  the  asymptotes  of  the 
original  hyperbolic  boundary,  into  which  this  boundary  deforms  as 
z-»0.  They  are  the  intersection  of  the  characteristic  cone 
through  the  point  xQ,yo  with  the  planforra.  The  Lorentz 
transformation  leaves  the  equations  for  the  hyperbola,  and 
therefore,  also  for  its  asymptotes  unchanged;  one  therefore  has 
for  the  other  boundaries 

Bn  +  5  -  xQ  for  n  >  0 

Bn  -  5  -  xQ  for  n  <  0. 

One,  therefore,  obtains  in  either  case 

;  .  Jo 

c  TpT 

or 
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*  n  p 

Now  the  expression  Qm  ,  Eq .  (62)  can  be  rewritten.  We  express 
dp/dn  by  Eq.  (65)  and  dFmn/dp  by  Eq.  (56).  Carrying  out  the 
integration  with  respect  to  p  last,  one  obtains 


Qmn2«  *a30“n  |1p[(a1-p)(a2+p)]~1/2fmn(p)Iin(p)dp 

“a2 

with 


(67) 


Im(p) 


(1+|P|) 

f 

"  I 


>2  ‘ 

,'2  ■'  2T2Tf72 


(1-0  cci-e)  -p  ri 


(68) 


One  notices  that  the  integral  im(p)  does  not  depend  upon  n.  The 
integrand  has  only  two  second  order  branch  points;  I ra ( p )  can 
therefore  be  evaluated  analytically.  The  integration  needed  in 
Eq.  (67)  is  best  carried  out  numerically. 

It  was  mentioned  above  the  values  of  m  are  determined  by  the 
character  of  the  upwash,  ra-O,  1,  2  will  probably  be  sufficient 
for  all  practical  purposes.  The  evaluation  of  the  integrals 
(68)  is  shown  in  Appendix  B.  Here  we  give  the  results  of  m-0,  1 
and  2. 


Let 


g(0V)  -  (1-P2r1/2 


g(1)(p2)  -  d-p2r3/2(3-2p2) 


(69) 


g(2)(p2)  -  (1-p2)_p/2((6-(15/2)p2  +  3p1<) 
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and 


g(0)  .0 


:(1)  -  -d-p2)-1 


i(2)  -  (-(7/2)  +  2p2)(1-p2)“2 


Im(p)  -  I®(p)  +  I®(P)  +  I®(P> 

im(p)  -  p~2  +  (p2) dogd  +(i-p2)1  /2)  +  i(ra) 


(70) 


l“(p>  -  -»  Sji  p'1  (71) 

l“(p>  -  -gm(p2)log|p| 

The  distinction  between  I®,  I®,  I®,  has  been  made,  because  they 

must  be  treated  separately  in  the  integration  with  re3pect^to  p. 

This  is  obvious  from  the  terms  p“2,  |p|~  and  log | p  |  in  Ia  , 

l(°),  and  I^,  respectively.  Accordingly,  we  write  for  the 
D  C 

expression  (67) 

~mri2  m  ran2  +  Qmn2  +  Qnn2  (72) 

w  b  c 

where  the  right  hand  sides  arise  by  replacing  I®(p)  by  the 
respective  expression  I®(p),  I®(p)  and  IC(P)*  ^ne  ^as 
specifically 


i3fTn  |Vai-P)(VP)]'1/2lP‘1+p^)log(^(H2)1/2)>Pg{m)]Ap)dp 


The  term  p_1  in  the  second  bracket  introduces  a  singularity  at  p-0. 
Further  singularities  are  encountered  at  the  upper  and  lower  limit 
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because  of  the  first  bracket.  The  occurrence  of  the  singularity  at 
p-0  is  to  be  expected  according  to  the  discussions  carried  out  in 
Section  VI  in  conjunction  with  Eq.  (25).  One  expects  that  at  this 
point  one  must  take  the  Cauchy  principal  value.  This  is 
demonstrated  in  Appendix  B.  Practically,  one  can  proceed  in  the 
following  manner. 

It  is  shown  in  Appendix  C  that 


1 

P  f  [(arp)(a2+p)]“1/2p"1dp»0 
~a2 

where  P  expresses  that  one  has  to  take  the  Cauchy  principal  value. 
Accordingly,  one  has 


,mn2 


■a38 


-n 


1 


[(a1-p)(a2+p)] 


'1/2J™(p)dp 


-a. 


(73) 


wi  th 

Tmn,  . 
J  (p) 


fmn(p)- 


fmn(0) 


+pg  ^  (p2)log(  1  +( 1-p2)1  /2)+pg^  (p2) 


(74) 


The  expression  — — ----- — is  regular  at  p»0;  fmn(p),  fmn(0)  and 
the  limiting  value  for  p-0  of  this  expression  (if  it  should  be 
needed)  are  found  Jn  Eqs.  (57),  (58)  and  (59).  The  expression 
Jmn(p)  is  regular  throughout  the  region  of  integration  -a2<p<a1 . 
Because  of  the  square  roots  in  the  bracket  of  Eq.  (73)  we  set 


P 


~  ^2  +  ^2 
"  P  2  +  ~— 


One  then  obtains 


(75) 
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(76) 


+  1 

nmn2  0-n  f , .  *2.-1 /2  Tmn, 

Qa  -  -a^B  j(1-p  )  J  (p)dp 

-1 

In  Jmn(p),  one  must  express  p  in  terras  of  p. 
Finally,  we  set 


p  -  sin  u 


Then 


-T 12 

-a38''n  f  J:nn(p)du 
-ir/2 


(77) 


(78) 


Here  p  must  be  regarded  as  a  function  of  u.  The  integrand  is  a 
periodic  function  of  u.  In  this  case  the  trapezoidal  rule  (with 
equal  intervals  in  u)  gives  very  good  results. 

In  I™(p).  one  must  take  the  occurrence  of  the  absolute  value 
of  p  into  account. 


1 

Q^n2  -  a3&  nir(  (ra+1  )/2)  j  sign  p  [(a^p)  (a2+p)  ]  1/2fmn(p)dp 

“a2 

0 

-a3e"nir((m+1)/2)[-J[(arp)(a2+p)]'1/2fran(p)dp 


+|[(a1-p) (a2+p) ]_1 /2fmn(p)dp 
o 


To  remove  the  singularities  for  p  -  a.  ,  and  p  -  -a,  one  might  apply 

mn  2  ^ 

the  transformation  used  in  Q™  again.  But  because  the  regions 

Ci 


-a2  <  p  <  0,  and  0  <  p  <  are  treated  separately,  it  is  probably 
convenient  to  make  separate  transformations  in  the  two  regions. 

Accordingly,  we  set  for  the  region  -a2  <  p  <  0 

p  *  p  +  Sg  »  p  *  p  -  §2  (79) 

"  <) 

and  subsequently  p  -  u  (80) 

and  for  the  region  0  <  p  <  a1 


p  ■  s,  ■  p  ,  P  ■  s1-  p 


P  -  u  ‘ 


Then  one  obtains 


V2 

a2 


-a36_nTr((m+1  )/2)[-J  (a^S^u2)"1 /2fmn(u2-a2)du 


-  - 


i  ,  2 » - 1 /2 „mn / ~  2,  .  -i 

+  j  (a^+a^u  )  f  (S.J-U  )duj 


Within  the  regions  of  integration  the  integrands  are  regular 
functions  of  u. 

The  expression 


Q®n2  -  -a3ffn  {[(a1-p)(a2+p)]"1/2p  g (m)  (p2)fmn(p)log  |  p  |dp 


requires  special  treatment  because  of  the  factor  log|p|,  although 
this  factor  is  not  detrimental  to  the  convergence  of  the  integral. 


A  small  region  -b  <  p  <  +b  around  p-0  is,  therefore,  treated 
separately. 

We  write 


q (p)  -  [(arp)(a2+p)]“1/2p  g(m)(p2)fmn(p)  (84) 


Then 


>1 


-b  +b 

3^n2-  -a3p"n[J  q(p )log | p |dp+ jq (p )log | p |dp+|  q(p)log|p|dp] 
~  ->b  b 

“3a 


(85) 


In  the  region  -b  <  p  <  b,  the  function  q(p)  is  developed  with 
respect  to  p.  Because  of  the  factor  p  contained  in  q  Eq.  (84), 
there  is  no  constant  term  in  the  development: 

2  3 

q (p)  -  t^p  +  &2p  +  b3  P  +*“ 


Consider  the  integral 


+b  +b 

f  (23 

j q (p )log  j p  |dp  -  j  (b,jp+b2p  +b^p  +  . . .  )log  |p  |dp 

-b  -b 

The  terms  with  odd  powers  of  p  vanish  because  of  the  antisymmetry 
of  the  integrand.  Taking  b  small  enough,  one  can  restrict  oneself 
to  terms  up  to  p  ;  then  only  the  coefficient  b2  is  needed.  One 
finds 


b2  -  gffl(0)[a1a2]“,,"[ 


1  /2rdf 


mn 


dp 


p-0 


1 ,1 
+  ?(~ 


‘1 


i_ 

a. 


)fmn(o)] 


(86) 
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are  found  in  Eqs.  (58)  and 


The  expressions  fmn(0) 
(59).  Then 


and 


+b  +b 

| q ( P ) 1 og | p | d p  -  b2  |  p2logJp|dp  -  ( 2b2/3)b3[(log  b)-(1/3) 3  (87) 

-b  -b 

To  remove  the  singularity  at  p  -  -a2  in  the  integral  for  the  region 
-a2  <  p  <  b  we  set  as  above 

p  +  a2  «  u2  ,  p  =  u2  -  a2  (88) 


Then 

(§2~b ) 1 /2 

2  |  (a1+a2-u2)“1/2fmn(p)g(ra)(p2)p  log  (a2-u2)du 
o 

(89) 

Here  p  is  regarded  as  a  function  of  u.  The  integrand  is  regular 
within  the  region  of  integration.  One  may  be  concerned  of  the 

fact,  that  b  is  small  and  that  in  the  vicinity  of  the  upper  limit 

“1/2  p 

u  «  (a2~b)  ,  the  argument  of  leg  (a2~u  )  is  small,  and  that, 

therefore,  |log(a2~u2)j  increases  rapidly.  This  would  make  small 

steps  in  the  numerical  integration  necessary.  A  smoother  procedure 

is  obtained  by  a  further  transformation 

log(a21 /2  -  u)  -  -v 

u  -  aj1 /2  -  exp(-v)  (90) 

'* 

du  -  exp  (-v)du 


-b 

|  q (p )log | p |dp  - 


H2 


1  /2 

The  limit  of  integration  u-0  becomes  v  -  -log(S2  )•  The  of 

integration  u  -  (S2-b)1/2  becomes  v  -  "log^S^^  ^Og-b)1  ^2].  Then 
the  integral  Eq.  (89)  is  transformed  into 

-log(S21/2-(S2~b)1/2 

2  j  (S1  +  a2-u2)fmn(p)g(m)(p2)  p  [log(a2/2+u]?1v]exp(-v)dv  (91) 

-log(S2/2) 

Here  p  and  u  are  regarded  as  functions  of  v.  The  region  of 
integration  would  tend  to  infinity  if  one  allows  b  to  tend  to  zero 
The  introduction  of  v  makes  it  possible  to  carry  out  the 
integration  with  a  uniform  interval. 

The  same  procedure  is  applied  to  the  region  b  <  p  <  §1 .  One 

sets 


-  p  -  u‘ 


a1  -  u 


(92) 


This  gives 


<art»>'2 

2 1  (a,*a2-u2)“1/2finn(p)  g(ol)(p2>p  log(a,-u2)du 


The  further  transformation 


.1  /  2 


log  (S1  -u)  -  -v  ,  u 


a] /2  -  exp(-v) 


(93) 


leads  to 


-Iok  (3  1/2-b) 

I  ,,,  .2,1  /2„mn,  ,  (m) 2 


1  /2 


I  (S^gj-u  )  f  (p)g  (p  )p[log(a1  +u]“V ]exp (-v )dv  (94) 


-iog(a)1 /2 


Again  u  and  p  are  regarded  as  functions  of  v. 

The  numerical  work  proceeds  as  follows:  the  shape  of  the 
airfoil  enters  the  problem  through  the  value  of  a^  and  a2  and 
through  the  given  upwash  i.e.,  through  the  value  of  m  and  the 


function  w(y/x),  which  usually  has  the  form  (Eq.  34)  .  Moreover, 

2  1/2 

the  Mach  number  M  and  with  it  0  -  (M  -1)  are  known. 

It  is  assumed  that  the  function  routines)  have  been  written  for 


fran(p)  , 

Eq.  (57) 

fmn(0)  , 

Eq.  (58) 

(dfmn |dp) |  , 

0 

Eq.  (59) 

g^m)  ,  and  g^ 

Eq.  (69) 

Of  course,  the  latter  expressions  need  to  be  programmed  only 
for  the  values  of  m  in  which  one  is  interested. 

In  a  collocation  method  one  chooses  a  number  of  values  of  c^. 
The  number  must  at  least  be  equal  to  the  number  of  values  n. 
Initially,  one  will  probably  cover  the  whole  range  of  possible 
values  of  ck,  -a2  <  ck  <  a1 ,  in  order  to  demonstrate  how  well  the 
upwash  condition  can  be  satisfied  by  the  use  of  a  restricted  number 
of  values  n. 

a 

For  these  values  of  ck  one  evaluates  wm(ck).  These  are  the 
inhomogereous  terms  in  the  system  of  equation?  obtained  from 
Eq.  (63). 

On  the  right  3ide  a  finite  number  of  values  of  n  is  admitted. 
In  their  choice  one  will  take  the  symmetry  properties  of  the 
airfoil  into  account;  for  a  symmetric  planform  and  a  symmetric 
upwash  one  will  have  even  values  of  n. 

The  further  description  refers  to  a  specific  value  of  ck.  For 
this  value  one  determines  c,  Eq .  (66),  a.  and  Eq.  (50),  and  a 0, 

~  i  •  £  3 

Eq.  (54).  On  this  basis  one  finds  Qm  ,  Eq.  (64).  The  quantity 

O 

Q  is  defined  in  Eq.  (74).  For  the  pivotal  values  of  u  needed 

a  a 

for  the  numerical  integration  one  finds  p,  Eq.  (77),  P  Eq. 

(75)  and  Tmn(p),  Eq.  (74).  The  quantity  Q^n2  iz  defined  in  Eq. 
(83).  With  fmn  available  as  a  function  routine  the  integrands  can 
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be  evaluated  immediately.  It  is  likely  that  one  will  use  different 
pivotal  points  for  evaluating  the  two  integrals. 

Qmn2  eXpPessecj  in  Eq .  (85)  as  the  sum  of  three  integrals. 

C * 

The  second  integral  is  evaluated  analytically  for  some  choice  of  b. 
With  b2  defined  in  Eq.  (86),  the  result  is  given  by  Eq.  (87).  The 
two  other  integrals  are  defined  in  Eqs.  (91)  and  (94).  One  chooses 
pivotal  values  of  v  and  evaluates  in  the  first  case  u  from  Eq .  (90) 
and  p  from  Eq.  (88).  In  the  second  case  u  is  defined  in  terms  of  v 
by  Eq.  (93)  and  p  by  Eq.  (92). 

One  then  has  Eq.  (72) 

Qmn2  „  Qmn  +  Q“n  ♦  Qmn 

w  wb  c 

,  "mn  ~mn1  A  *  mn2 

and  Q^c  )  *  Q  +  Q 

k 

a 

The  Qmn(ck)  are  the  matrix  elements  in  the  system  arising  from 
Eq.  (63),  n  and  k  are  respectively  the  row  and  column  indices. 

After  the  coefficients  bn  are  found  from  this  system  of 
equations  one  can  determine  h(x,y)  from  Eqs.  (39)  and  (40). 
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SECTION  VIII 


IMPLEMENTATION 

In  many  cases  the  planform  has  a  tip  in  which  two  subsonic 
edges  meet.  Then  part  of  the  flow  can  be  computed  as  a  conical 
field,  for  instance  by  the  procedure  shown  in  Section  VII.  One 
determines  the  coefficients  bn.  Then  Eq.  (35)  gives  an 
analytical  approximation  for  the  flow  field.  From  these 
expressions  one  can  compute  a  pointwise  representation  of  h. 

Beyond  the  region  for  which  the  evaluation  as  a  conical 
field  is  possible  one  uses  Eq .  (32)  to  determine  h  at  prechosen 
points  of  the  planform,  which  will  be  called  control  points.  We 
assume  that  these  points  lie  on  lines  £  -  const.  Since  the 
boundaries  of  the  regions  of  integrations  are  Mach  lines,  it 
seems  reasonable  to  impose  the  further  condition  that  the  points 
lie  on  preselected  Mach  lines.  One  thus  obtains  an  arrangement 
similar  to  Mach  boxes,  but  without  the  assumption  that  within  a 
Mach  box  the  potential  is  approximated  by  a  constant.  The  author 
rather  thinks  in  terms  of  a  pointwise  representation  of  the  flow 
field.  The  function  h  at  a  certain  point  (x,y)  Is  determined 
only  by  data  at  the  planform  within  its  forecone.  The 
computation  can  therefore  be  carried  out  by  marching  in  the  x 
direction.  In  evaluating  the  integrals  in  Eq.  (32)  it  is  assumed 
that  the  integrations  with  respect  to  n  are  carried  out  first. 
Because  the  computation  marches  in  the  x-direction,  the 
integrands  in  the  integration  with  respect  to  n  are  known,  except 
for  the  immediate  vicinity  of  the  point  (x,y)  under 
consideration.  From  the  values  of  h  at  the  control  points  one 
can  derive  an  approximation  for  intermediate  points.  One  r<ght 
consider  an  approximation  that  is  piecewise  linear  in  n,  but  it 
is  probably  preferable  to  use  a  higher  order  interpolation 
because  of  the  denominator 

[(x-£)  -$  (y-n)  1  which  becomes  zero  for  (x-£)  =  +_  B ( n~ y ) . 

One  has  the  choice  of  using  information  for  h  that  comes  only 
from  within  the  region  of  integration,  or  one  may  involve  also 
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points  from  outside  if  one  can  count  on  the  smoothness  of  h.  The 
first  possibility  may  be  simpler. 

We  have  seen,  that  in  the  vicinity  of  the  point  (x,y)  the 

integral  over  n  does  not  vanish,  although  at  the  point  (x,y)  the 

upper  and  lower  limit  of  the  integral  coincide.  A  limiting  value 
(22) 

h  -  h  is  encountered  for  n  *  y.  But  at  the  station  (x,y) 
and  the  station  upstream  of  it  this  limiting  value  is  not 
available  (unless  one  admits  information  from  outside  the  region 
of  influence).  Here  the  following  procedure  is  suggested.  We 
assume  that  in  the  entire  triangle  shown  in  Figure  5,  h  is  given 
by  one  analytical  expression. 

The  starting  point  is  Eq.  (31).  Developing  the  integrand 
with  respect  to  n  one  obtains 


„-if  rn2ht22)(S,y)detln  „-lL22 

*  JJ-TTW7*  8  Jh 


f  "2  - 

I  n  dn 


h~U,y)(  ) 

J  (i-n  )'  /2 


B1  (ir/2)  fh{22)  U,y)d£ 


Along  the  line  2,0,1  one  has 


h.  +h0-2hn 
1  2  o 


where  An  is  the  interval  in  n-  In  the  evaluation  of  the  flow 
field  one  obtains  in  an  intermediate  step  h^.  Therefore,  h^  is 
available  at  points  2,0,1  and  one  has  also 

t.u.h25-2h05 


For  the  extrapolation  one  therefore  can  use 


h1+h0-2ho  h.  c  +h.,_-2h,_ 

h22(y,5)  -  -1— - (5-0  -^-^-2— 

*  An  An 


Carrying  out  the  integration  from  C  «  £0  to  C  -  5S  *  £0  +  BAn  one 
then  obtains  for  the  contribution  of  the  triangular  region  to  the 
integral 


tt/2  [(h1+h2-2h0)/An)  +  (8/2)(hu+h2c-2ho5)] 


where  6n  is  the  interval  in  n. 
Furthermore 


nn5 


h£1 +h52“hco 

V 


This  is  used  for  an  extrapolation 


h,  +h5-2h  h,.,  +h^^-2, 

h(2‘U>  -  —  -  ♦ 


An 


T‘ 


nr1 +nco-dCn 


An 


The  interval  of  integration  extends  from  £<»0  (at  points  0,12)  to 
£  «  8An  at  point  5. 

The  contribution  of  the  entire  triangle  is  therefore  given 
by 


11/2  h  [h1+h2”2ho  +  I  Ari(hu+h2r2V)] 
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APPENDIX  A 


DISCUSSION  OF  DERIVATIVES  OF  <|>(3)  WITH  RESPECT  TO  z  WITHIN  AREA  II 

The  area  II  is  defined  in  the  main  text.  It  is  bounded  by 

the  leading  edge,  by  the  straight  line  5  (where  ^  is  chosen 

so  that  this  line  intersects  the  hyperbola  (1-?)2-(n-y)2  -  z2  =  0 

twice  within  the  planform)  and  by  two  portions  of  this  hyperbola. 

One  notices  that  the  variable  z  with  respect  to  which  the 

s  T I 

differentiations  are  carried  out  enters  the  expression  <{>  only 

2 

in  the  form  of  the  parameter  s  =  z  . 

We  write 


(x ,y ,z ) 


4>(x,y  ,s) 


Then 

4>(dII)  (x.y  ,z)  =  2z<(>s(x  ,y  ,s) 

<J>2dII)  (x.y.z)  =  20s(x,y,s)  +  4z2(j)gg(x  ,y  ,s) 
In  the  limit  z  =  0  one  obtains 

^(d,11) (x  >y >z)  „  o 

JI)  (x  ,y  ,z)  -  <|>3(x,y,s) 


provided  that  <b  exists.  Whether  *  exists  in  the  limit  z  -  0 

3  8  S  S 

is  not  immediately  obvious,  this  is  the  reason  for  the 
discussions  carried  out  in  this  Appendix.  But  after  the 
existence  of  <}>gs  has  been  established,  it  is  no  longer  necessary 
to  evaluate  it  in  detail,  because  we  are  only  interested  in  the 
limit  z  -*>  0. 

The  leading  edge  is  given  by  the  equation 


n 


O  V  t 


( A .  1  ) 
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The  inverse  of  g  is  denoted  by  f 


C  -  f<n) 


(A. 2) 


We  have  introduced 

-  r,  _ >  2  .1/2 

n  =  y  +  n  [(x-O  -s] 

n  =  ( n-y ) [ (x- c)2-s]  1  /2 

In  the  Cn-plane  the  leading  edge  is  then  given  by 

n  =  g  ( £ ,  s  ,x  ,y )  (A.  3) 


with 


gU.s.x.y)  =  (g(£)-y)[(x-£)2-s]*1/2  (A.  4) 

In  the  present  context  the  arguments  x  and  y  are  kept 
constant.  In  many  cases,  they  will  therefore  not  be  listed  as 
arguments  in  the  functions  in  question.  The  inverse  of  the 
function  g  at  constant  s,  x,  and  y  is  denoted  by  f(n,s,x,y). 
Accordingly,  we  have  as  alternative  to  Eq .  (A. 3)  for  the  equation 
of  the  leading  edge  in  the  £n-system 

5  =  f(n,s,x,y)  (A. 5) 

For  subsonic  leading  edges  one  has  |dg/d£|  <  g”1.  For 

simplicity,  the  present  discussions  will  be  carried  out  for 
2 

M  =  2.  Then  8=1.  Depending  upon  the  orientation  of  the 
portion  of  the  leading  edge  under  consideration  the  function  g  is 
either  monotonically  increasing  or  decreasing.  For  g,  however, 
this  is  not  necessarily  correct.  In  Figure  6,  the  curve  AB  is  a 
leading  edge  along  which  n  decreases  monotonically  with  £;  yet  it 
is  tangent  to  a  straight  line  through  the  point  (x,y).  Such  a 
straight  line  is  the  limit  of  some  curve  n  =  const  as  z  *  0.  The 
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variable  n  therefore  first  increases  then  decreases  as  one 
travels  from  A  to  B  along  the  leading  edge.  If  this  should  occur 
the  area  is  denoted  by  lib;  in  the  ^-directive  it  extends  from 
to 

Let  us  first  consider  the  problem  with  subsonic  leading 

edges.  Figure  7  shows  such  a  configuration  in  the  5n  -  and  in 

the  5n-planes.  At  the  leading  edge  the  potential  and  with  it  h 

have  a  square  root  singularity.  The  detailed  discussions  carried 

out  here,  are  made  necessary  because  the  leading  edge  moves  in 

the  5n  system  as  s  varies.  We  denote  by  y  the  local  sweep  angle 

and  introduce  a  local  Cartesian  coordinate  system  uv ,  where  the 

1  1 

u-direction  is  perpendicular  to  the  leading  edge.  Let  £  ,  n  be 

the  coordinates  of  the  point  at  the  leading  edge  under 

consideration.  Then  (Fig.  8) 

=  u  cost  +  v  sinY  u  =  (S-^^cosY  -  (n-n1)  sinY 

n-n1  -  -u  sinY  +  v  cosY  v  =  (S-i^sinY  +  (n  V)  cosY 

Because  of  the  square  root  singularity  one  has  locally 

h  -  const  u1/2  =  const[(  )cosY  -  (n-n1  JsinY]1  ^2 

Therefore  for  n  =  n1  =  const 

h  =  const(cosY) 1 /2(^-^1 ) 1 /2 

used  for  £  =  ^  =  const 

h  =  const  (sinY)1^2  (V-n)1^ 

The  flow  field  in  the  vicinity  of  the  intersection  of  the 
two  subsonic  leading  edges  is  rather  complicated.  The  basic 
structure  is  that  of  a  conical  field,  but  higher  order 
corrections  occur  if  the  leading  edges  are  not  straight  or  if  the 
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upwash  is  not  constant.  This  vicinity  is  best  considered  in  the 
original  £n-system.  This  is  the  region  Ha  in  Figure  9.  In 
order  to  avoid  dealing  at  the  same  time  with  two  boundaries  that 
depend  upon  s,  the  remainder  of  the  region  is  subdivided  by  a 
line  n  -  0.  One  thus  obtains  the  regions  lie  and  lid.  In  the 
vicinity  of  the  point  where  a  line  n  =  const  is  tangent  to  the 
leading  edge,  if  such  a  point  should  occur,  we  introduce  for  this 
vicinity  a  separate  region  lib. 

In  the  tip  region  Ila,  uhe  boundaries  do  not  depend  upon  s. 

One  has 


<t>(s,Ha)(x,y,z) 

<J>(d’IIa)(x,y,z) 


ff 

JJ 

zf 

ZJ 


hU,r))d£dri 


[(x-C)2-(ri-y)2-s]1  /2 

h(5il0dCdn  _ _ 

[(x-0?-(n-y)^s]3/2 


The  denominator  is  always  different  from  zero.  Hence 

<j>(d,IIa)(x,y,0)  =  o 


(d  ,  Ila) ,  v  n,  ff 
<p  U,y,0)  «  j  j 


hU.n)<Udn 


n-y)2]3/2 


( A .  6) 


To  bring  the  expression  into  a  form  which  will  be  encountered  in 

ether  regions,  we  introduce  n.  Specializing  immediately  to 
2 

s  =  z  =0  one  has  n  =  y  +  n(x-£). 


Then 


A  ( d  ,  1 1  a ) ,  „  ^  f  fh(C,y  +  h(x-5)  )d£dn 

U,y,0)  -  j  j  — - — r  j~T7T ' 


(x-C)  (1 


■n2>3/2 


Now 


-  (  .  -  9.  1  /?  )  =  T  -9 


dn  (l-n2)1 /2  <1-n2)3/2 
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In  the  last  equation  we  carry  out  the  integration  with  respect  to 
n  first,  and  transform  the  inner  integral  by  an  integration  by 
parts.  At  the  upper  and  lower  limits  (the  intersections  of  a 
line  £  =  const  with  the  leading  edges)  h  vanishes  and  one  obtains 


4d,IIa>‘x.y.»)  ■  -  (|1'(2.>.U.y^(«-c)..).ndn)ds  U-7) 


(2)  -  -1/2 
The  function  h  behaves  at  the  limits  as  An  ;  the  inner 

integral  converges.  The  denominator  does  not  vanish  since 

n  <  1.  Returning  to  the  original  coordinates  one  obtains  an 

integrand  familiar  from  Eq .  (32) 


<^d,IIa)(x,y,0) 


ff 

JJ 


(n-y)h(2)  U,n) 
- -2— - 


Alia  (*-5)  Ux-O  -(n-yri 


(A. 8) 


p 

One  remembers  that  in  this  Appendix  g  =  1.  In  a  region  lib 
(if  it  should  occur)  the  integration  with  respect  to  n  is  carried 
out  first.  We  write 


♦ 


(s) 


f 


j  FU,s)d£ 


(A. 9) 


where 


n=gU,s) 

FU.s)  =  [  dr, 

J0  ( 1-n) 1 

qU,x,s)  =  ((x-^)2-s)' /2 


s  =  z 


At  the  upper  limit  (which  is  a  point  of  the  leading  edge)  h  has  a 
square  root  singularity 
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h  -  const  (g(S)-n)1 /2 

This  expression  is  rewritten  in  terras  of  n  and  g,  Eqs.  (13)  and 
(17) 

—  I  /p 

h  -  const  [g(5)-y-n*q] 
h  -  const  q  [gU)~n)  l1  /2 

The  factor  q(5,s)  has  no  singularity  in  the  whole  of  the 

2  2 

region  II.  The  discussions  of  Appendix  E  show,  that  8  F  /9s  is 
bounded.  One  therefore  finds 


Here 


9F/3s  |  = 
s  *=  0 


(1/2) 


n-g(e.o) 

f  nh 


J 


(2)  (C , (y +n(x-g ))dn 

7x-o)r-?)172 


Returning  to  the  original  coordinates  £,n,  one  obtains  the 
integrand  occurring  in  Eq.  (32) 

The  regions  lie  (or  lid)  are  subdivided  into  regions  lid 
and  IIc2  by  the  line  n=g(52»s)»  Figure  11.  This  boundary  moves 
as  s  changes.  Consider 


4><Hc1>( 


x.y.z)  = 


F( n,s)dn 

-27T72 
_  _(  1-ti  ) 

n=iU2’s) 


(A. 10) 


where 


5>i 


(A. 11) 


^1 

F(n,s)  »  |  h(£,(y+nq(£,x,s))d£ 
f (n, s) 

As  before  f(n,s)  denotes  the  inverse  of  the  function  n»g(£ ,x ,y ,s) 
at  constant  x,y,s.  The  lower  limit  in  F  lies  at  a  point  of  the 
leading  edge.  There  the  function  h  has  a  square  root 
singularity.  We  write 

h(£,y+r,q)1/2  -  U-f ( n , s ) ] 1  /2h(£, n,s) 

(The  dependence  upon  x  and  y  does  not  appear  in  the  arguments 
because  they  are  kept  constant.)  Here  h  is  free  of 
singularities.  Now  the  discussions  of  Appendix  E  are  again 
applicable.  The  first  and  second  derivative  of  F  with  respect  to 
s  are  bounded,  and  as  the  second  derivative  wit"  respect  to  z  is 
needed  for  z  =  0.  Only,  there  is  no  need  to  e*0.uate  Fss' 
obtained  from  Eq.  ( A . 1 1 ) .  At  the  lower  limit  the  function  h 
vanishes.  The  discussions  of  Appendix  E  show,  that  (for  the 
first  derivative)  the  square  root  singularity  encountered  at  this 
point  does  not  matter. 


F 


s 


,,,  f  nhl2> (5.y*nQ(5,x,s))dn 

(,/2)  j - - - 

f (n,s) 


(A. 12) 


Then,  by  differentiating  Eq .  (A. 10)  with  respect  to  s  and 
substituting  Eq .  (A. 12). 
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^ (s , lib ) 
s (x ,y ,0) 


Og/3s) 

5 


I  F(g(c2,0),0) 
ljl-i(52.0)"2)W5i 


3*=0 


n 

-27X71 


( 1-n  ) 


n-gU2,o) 
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(j  h(2) U,y+n(x-Q)dg)d- 
f  (ri,0) 


(A. 13) 


Going  back  to  the  original  coordinates  5  and  n,  one  obtains  again 
an  expression  with  the  integrand  of  Eq.  (32),  except  for  the 
first  term  on  the  right,  which  arises  because  of  the 
s-dependence  of  the  lower  limit  in  Eq .  (A. 10) 


In  the  region  lie  one  has 


n=gU2,s) 

♦  (s'IIo)(x,y,z)  -  f  (A. 11) 

l  (1-n  ) 1 /2 

with 

51 

F(n,s)  *=  j  h(C,(y+nq(C,x,s))dS  (A. 15) 

^2 

*1 

=  |  U-f(n.s))1 /2hU,n,s)d5 
^2 

A  comparison  of  Eq.  (A. 15)  and  (A. 11)  shows,  that 

F(g(£2,s),s)  =  F(gU2,s),s)  (A. 16) 

The  limits  in  the  function  F  do  not  depend  upor  s.  A 
complication  arises,  however,  because  for  n=g(C2*s^  the-  inteSrand 
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has  a  point  in  common  with  the  leading  edge.  One  will  surmise 
that  this  is  unessential,  but  some  discussion  is  desirable.  One 
obtains  from  Eq .  (A. 15) 

Fs(n,s)=|(-1/2)[f(2)(nfs)(c-f(n,s))'1/2h(c,n,s) 

^2 

+  U-f  (n.s))1  /2h(3)  (5,n,s)]de  (A. IT) 

The  first  term  in  the  integrand  tends  to  infinity  for  the  lower 
limit  S-S2  ^  h-gCC^s),  that  is  at  the  upper  limit  in  Eq.  ( A . 1 3 ) 
but  the  integral  converges.  In  addition,  one  must  show,  that  F 

ss 

is  bounded  (see  Appendix  E) .  One  has 

F33(n,s)  -  |{[-(l/i|)(f(2)(n,3))2(c-f(n,s))“3/2 
^2 

-(i/2)f(2*2)(5-f(n,3))"1/2]h(5,n,s) 

-f(2)(n,s)U-f(n,s)~1/2h(3)U,n,s)  +  U-f(n,s)1/2h(33(£  ,n,s)}dl- 

After  the  integration  has  been  carried  out  the  term  with  the 
factor  (5~f(ns))  gives  a  result  which  tends  io  infinity  as  n 
approaches  iU2,s).  This  contribution  is 

0(c2-f(n,s))"1/2 


Let  n=gU2,s)-e  (A. 18) 

Then  f(n,s)  =  f (gU2,s)-e ,s) 

But  f  is  the  inverse  of  g,  therefore 
f  (gU2,s)  ,s)  =  i2 
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and 


f (g2,s)-e .s)  “  S2  ~  °( e  ) 


The  critical  tern 


U2~?  <*>, 

(  0  T  T  /->} 

If  one  forms 

3  3 

of  n<»gU2,s) . 


f 

j 


0(  e 


-1/2 


- 1  /?  - 1 /? 

s )  is  therefore  0^ e 

from  Eq.  (A. 13)  one  obtains 
)de  -  CKe1^2) 


). 

in  the  vicinity 


With  e  defined  by  Eq .  (A.18)  the  contribution  of  F  to  a  is 

S3  3  S 

therefore  bounded.  It  suffices,  therefore,  if  one  evaluates  only 


. (s,nc) 


Og/3s)|(1-£U2,0)2)"1/2F(£U2,0)  ,0) 


3-0 


+ 


n-gU2 

J 

0 


^  (Ti,0)dn 

(1-n2)1/2 


(A. 19) 


Because  of  Eq.  (A.16)  the  first  term  on  the  right  in  this 

expression  cancels  the  first  term  on  the  right  in  Eq.  (A. 13)*  In 

~( 2) 

the  second  term  in  Eq.  (A. 19),  F  =  F  is  substituted  from  Eq. 

s 

(A. 17),  furthermore,  h  (instead  of  R)  is  reintroduced. 

Returning  to  the  coordinates  £  and  n  one  arrives  again  at  the 
integrand  in  Eq.  (3?). 

If  the  leading  edge  is  partially  subsonic  and  partially 
supersonic  one  can  distinguish  on  the  wing  surface  between  the 
regions  A  and  B  (Fig.  12).  They  are  separated  by  the  Mach  wave 
emanating  from  the  corner  between  the  subsonic  and  supersonic 
leading  portions  of  the  leading  edge.  The  boundary  of  the  region 
of  integration  for  a  given  point  (x,y,z)  always  consists  of  the 
hyperbola  and  part  of  the  leading  edge.  If  the  point  (x,y)  lies 
in  the  region  A,  then  the  part  formed  by  the  leading  edge  is 
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entirely  supersonic.  At  the  leading  edge  the  potential  and  with 
it  h  are  zero,  as  always,  and  h  increases  linearly  with  the 
distance  from  the  leading  edge.  The  precautions  which  we  took 
because  of  the  square  root  singularity  at  a  subsonic  leading  edge 
and  because  of  the  conical  field  at  a  tip  formed  by  two  subsonic 
leading  edges  are  not  needed,  and  one  arrives  without  the 
detailed  discussions  shown  for  such  cases  in  this  appendix  at 
Eq .  (32).  If  no  points  other  than  those  of  region  A  were 
considered,  it  would  be  preferable  to  represent  the  flow  field  by 
a  source  distribution.  The  function  h  is  then  directly 
determined  by  the  given  upwash  Eq.  (23)  and  the  potential  is 
found  by  integrations,  for  instance  from  Eq.  (18).  The  flow 
field  ir  the  region  A  is  not  affected  by  the  conical  field  which 
arises  at  the  juncture  of  the  supersonic  and  subsonic  leading 
edges . 

The  cases  where  the  point  (x,y)  lies  within  the  region  B 
differs  from  those  with  two  subsonic  edges  by  the  treatment  of 
the  vicinity  of  the  juncture  of  the  two  edges.  Figure  13.  We 
distinguish  between  the  subregion  CDE,  lying  downstream  of  the 
Mach  wave  CD  (Region  Hal),  and  the  subregion  CDFG,  downstream  of 
the  supersonic  part  of  the  leading  edge,  but  upstream  of  the  Mach 
wave  CD  (Region  Ila?) . 

The  character  of  the  flow  field  in  the  region  Hal  is 
similar  to  that  of  the  region  formerly  denoted  by  Ila.  The 
contribution  of  the  region  Hal  to  the  upwash  at  the  point  (x,y) 
is  again  evaluated  in  the  original  Subsystem,  and  one  obtains  as 
an  intermediate  result  Eq.  (A. 6).  In  a  further  step  performed  to 
bring  this  result  into  the  form  which  arises  in  other  regions,  we 
introduced  n  instead  of  n  and  carried  out  an  integration  by  parts 
Eq.  (A. 7).  Subsequently,  we  returned  to  the  original  £n  system 
Eq.  (A, 8).  This  procedure  is  applied  again.  But  here  the  lower 
limit  for  n  is  not  a  leading  edge  where  h  is  zero,  but  the  Mach 
wave  CD,  where  the  potential  is  different  from  zero. 
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One  obtains 


4> 


(dlla) 


[? _ h(s.n)  (n-y) _ 

1  ( x- e )  (x- 5 ) "2-  ( n- y ) 2  ] 1  ^ 


rJl 


d£ 


-  ffr 
iJ[' 


h(2)  (S,n)  (n-y)d£dn 


n-n-U- 


1 1 A 1 


(x-C)  C (X 


C)2-  (n-y)2]1  /2 


(A20) 


The  first  term  on  the  right  is  the  contribution  of  the  Mach  wave 
CD.  The  second  terra  has  the  form  familiar  from  previous 
discussions.  In  the  first  term  we  introduce  n  instead  of  £  as 
variable  of  integration 


..  n  **  Dq  -  (6-6q) 

One  remembers,  that  for  Mach  waves  d£/dn*+1 ,  because  here  B«1 . 
Then  the  integral  assumes  the  form 


n 

r 

J 


c 

( 


nc 


hU,n)  (n-y) 

2'- 


(x-02C(x-02-(n-y)2]1/2 


>1 


dn 


S=6C-  ( n- nc) 


(A. 21 ) 


We  shall  show  that  this  term  will  be  cancelled,  provided  that  one 
expresses  the  solution  in  the  region  IIa2  by  a  doublet 
distribution,  and  evaluates  the  upwash  at  the  point  (x,y) 
accordingly.  It  is  true,  the  determination  of  the  solution  in 
region  A  is  simpler  if  one  applies  a  source  distribution,  the  use 
of  the  doublets  has  the  advantage  that  one  obtains  in  all  case 
the  same  expressions  in  the  integral  equation. 

Because  part  of  the  boundary  of  the  region  IIa2  is  formed  by 
the  hyperbola  we  introduce  again  n  and  perform  the  integration 
with  respect  to  £  first.  Then 
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♦ 


(s,IIa2) 


Tip 

f  F(n.s) 

- - ? TT>  dn 

J  ( 1-n2) 1 /2 


(A. 22) 


n--i 


wi  th 


5upper(^ 

F  (  ri ,  s  )  -  j  hU,y+nq)di- 
f (n,s) 


Here,  since  the  discussions  are  carried  out  for  g-1 


q  -  C(x-5)2-s]1 /2  ,  s«z2. 


For  h-Hq  the  two  limits  of  the  last  integral  coincide. 
Therefore , 


F(nc,s)  =0 

The  function  E>f(n,s)  is  the  equation  of  the  (supersonic)  leading 
edge  in  the  system. 

The  upper  limit  6upper  is  partially  given  by  the  line  DF, 
and  partially  by  the  line  CD. 

Differentiating  Eq .  (A. 22)  with  respect  to  z,  one  obtains 


Now 


$ 


(d , IIa2) 


f  F(2)(n,s)d£ 


'  J 

n=-l 


-2x1/2 

(1-n  ) 


F(2)(n,s) 


^hU,y*nq)| 


upper 


(A. 23) 


(A. 24) 
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The  last  term,  evaluated  for  s»0  and  substituted  into  Eq.  (A. 23) 
gives  a  contribution 


f  fnhv2^  U,y+n(x-0 

ii  d--,27^ir 


d^dn 


In  terms  of  the  original  coordinates  5n  it  assumes  the  form 
familiar  from  Eq.  (32). 

In  the  first  term  of  Eq.  (A. 24)  dSupper/ds=.0  along  DF. 

Since  the  slope  of  CD  for  3=1  is  -1,  one  has  in  the  £n  plane 


n  =  y  +  nq 

Then  Eq.  (A. 25)  for  5upper  becomes,  after  substitution  of  q 

*  (x-«c>  *  (*-*o>  *  "[<x-«Upper)2-s]1  <A-26) 

We  note  that  one  obtains  for  s=0 

■  (x"5upper)  (1"^  +  (x'xc)  +  (y_yc)  =  0  (A. 27) 

To  obtain  d£upper/ds  we  differentiate  Eq .  (A. 26)  at  constant  ii 
with  respect  to  s 
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-(x-C  )dE  - ( 1 /2)ds 
ric  +  „  __ _ upper'  ^upper 

UPPer 


and  after  specialization  to  s=0 


dCnnno  /ds 
upper 


1 

2 


(A. 28) 


<x-5upper)<1-,') 


This  is  now  substituted  into  the  first  term  of'  Eq .  (A. 24)  and 
subsequently  into  Eq.  (A. 23). 

Then  one  obtains 


ch u 

[  (x-c; 


f  n^upper(^)>y" 

( n) )  ( 1  -n) 


'D 


upper 


To  compare  this  expression  it  with  Eq.  (A. 21),  one  must  replace 
the  variable  n  by  n.  One  has  from  Eq.  (A. 14.3) 


d5upper(1  n)  *  <x-«upper)dr|-0 


Along  CD,  <Hupper  =  -dn.  Therefore 


dn  =  ( 1-rj) 

dn  x-5 


upper 


Thus  the  above  expression  assumes  the  form 


hU,n)  (n-y) 


(x-02C(x-02-(n-y)2]1/2 


dn 


^^c”(n"ric) 


Because  the  limits  of  integration  are  interchanged,  this  is, 
indeed,  the  negative  of  the  expression  (A. 21).  In  the  combined 
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contributions  of  the  regions  Hal  and  IIa2,  the  integrals  along 
the  line  CD  cancel . 


APPENDIX  B 


EVALUATION  OF  CERTAIN  INTEGRALS 
The  integral  Im(p)  is  introduced  in  Eq.  (68) 


.  >.2‘ 

I  (P)  -  • - ^ - ^ 


i(l 


(i-c)2C(i - e) 2-P2i D1  /2 


( B.  1 ) 


The  radicand  is  rewritten 


(l-O2  -  P 2l2  =  ( 1-p2)-1  [( C(  1-P2)-1  )2-p23 


We  introduce 


C(1-p  )-1  -  -q 


(B.  2) 


Then 


A  4 

£  - 

1-P 


1  -  5 


SlPL 


dC 


1-p‘ 
~dq 


1"P‘ 


The  lower  limit  of  Im(p)  transforms  into  q=1  ,  the  upper  limit 
into  q  -  | p | .  Then  one  obtains 


Tm  (  \ 
I  (p) 


(1-P2) 


-m-(1/2  f  (1-q)m*2dq 
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One  has 


(4  \  ni+  2  r/i  2  \  f  2\nfn+2 

(1-q)  *  [( 1-p  )  -  (q-p  ) J 

, «  2sm+2  ,m+2w,  2,m+1  ,  2.  ^  >ra+2Wl  2Nm-2/„  „2s2 

-(1-p  )  -  (  1  )(1-p  )  (q-p  )  +  (  2  ) ( 1 ~P  )  (q-p  ) 


In  particular  for  m=2 


(1-Q)2  (1-P2)2  2( 1-p2)  „  . 

' 


for  m-1 


(1-q) 

(q-p2)' 


-ill #  -  *  (3-2p2)  -q 

(q-p2)2  (q-p2) 


for  m-2 


*  (6-8p2*3p4)  *  (-*1+2p2)q  .  q2 

(q-p2)  (q-p  )  (q-p  ) 


These  values  of  ra  are  probably  sufficient  for  practical  purposes; 
expression  for  greater  values  of  m  are  easily  derived  when 
needed.  We  introduce  the  following  integrals 


1-2  \ 


, J  ,  ,  2,2,  2  2.1/2 

I p |  (q-p  )  (q  -p  ) 


:-if 


i  .  ,  2s  ,  2  2x1/2 

| P |  (q-p  ) (q  ~P  ) 


:  -  I  =  f 

o  o  J 


(q2-p2)1/2 
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n 


|P 


qndq _ 

(qV)1'2 


Then 

I°(p)  -  (1-p2)3/2  I_2  -  2(1-p2)1/2  I_1  +  (1-p2)“1/2Io  (B.3) 

^(p)  -  (1-p2)3/2  I_2  -  3 ( 1  - p2) 1  /2  I_1  +  (1-p2)“3/2(3-2p2)Io 

-  ( 1 -p2 ) ~ 3 /21 1 

I2(p)  -  (1-p2)3/2  I _2  «  4(1-p2)1/2  I_ 1  ♦  (1-p2)“5/2(6-8p2+3p4)I0 

+  (1-p2)"5/2(-4+2p2)I1  +  (1-p2)“5/2I2 

The  integrals  In  (in  their  indefinite  form)  are  connected  by  the 
recurrence  relation 

In  -  (1/n)[qn'1(q2-p2)1/2  +  <n-1)p2In_2]  n  >  2 


This  is  shown  as  follows 


n 


r  n . 

I  Q  dq 


[  q"-1  q<iq 

J(q2-P2)’/2 

qn_1 (q2-p2)' /2-(n-1 ){qn‘2(q2-P2)’ /2dq 


,n-1(,V)1/2-(n-1)f 


i  qndq _  n.ti  q  U4 

'(q^W7*  (  PJ(q2-p2)’/2 


2  f  qn2dq 


The  two  integrals  on  the  right  are  respectively  In  and  In_2* 
This  leads  to  the  above  recurrence  relation.  One  verifies  by 
differentiation 
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I.  -  i . *  =  log  [q  +  (q2-p2)1/2]  -  log 


f  dq 

0  “  J  (q2-p2)T/2 

|P| 

1  1 
t  f  qdq  t  „ 2  _2J  /2  | 

1  '  i('^>  ‘  u  'p  5 


|P| 


IPS 


(1-P2)1/2 


u(i-p2)1/2 

Tp T 


IpI 


Then,  according  to  the  above  recurrence  relation 


I2  *  (1/2)  (1-p2)1/2  +  (P2/2)Iq 


The  integrals  I_2  and  I_^  can  be  found  in  tables.  The  author  has 
used  Ref.  1,  Formula  234. 3a  for  I_2  and,  with  some  modification, 
formula  234.3b  for  I  , .  One  has 


,  _  2  2n-1/2  _  a2-aq 

2TZ2^J72  -  (a  5  arc  003  iTF^y 


dq 


J (q-a) (q  -a  ) 
f  dq  ,  2  2.-1  [-(q2-a2)1  /2  .  f 

jT— -2-2-T7tT2  ■  U  -«  5  VS - -  “J 


(q-a)  (q  -a  ) 

At  present  a=p,  a^p*1,  then  a^-a^  «  p^(l-p^) 


«,-a)(q2-a2)'/2 


Therefore 


1 

I_1  =  |  p  j  _1  (1-p2) 1  /2  arc  cos 

Q_P  q-|P| 


IpI'1  (i-P2)"1/2Tr/2 


2  2  1  /2  1 

i.9-  p"2o-p2)-'  I  ♦  p2  i‘. 


q-p 


( B.  4) 


<H?I 


*-?  =  ~2  ( 1  _  p2  )  3  /2  t1+|p|ir/2] 

Substituting  these  expressions  into  Eq.  (B.3),  one  obtains 
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l0^=p^  ~  W  5  *  (1-P2)  1/2  [l0g(1+(l-p  ) 1  /2 )  -  log  I  p  I  3 

ll(p)— j  -  IT  +  ( 1-p2)~3/2( 3-2p2) [log( 1 +( 1-p2)1  /2)-log |p | ]-( 1-p2)-1 

I2(p  )  — 1  -  +  ( 1-p2)-2  (- (7/2)  +2p2 )  + 

+(1-p2)'5/2(6-(15/2)p2+3p4)Clog(1+(1-p2)1/2)-log|p| ] 


69 


APPENDIX  C 


TREATMENT  OF  A  CERTAIN  INTEGRAL 


The  integral 
31 

(i/p)[(arp)(a2+p)]“1/2fmn(p)  dp 

*2 

has  a  singular  integrand  at  p»0.  In  the  main  text  it  is  stated 
that  one  should  form  the  Cauchy  principal  value.  Here  we  shall 
justify  this  procedure.  The  term  arises  in  the  evaluation  of 
»  Eq.  (67);  specifically  it  is  due  to  the  function  Im(p)  , 
defined  in  Eqs.  (68)  (70)  and  (71).  The  above  expression  arises 
from  the  first  term  on  the  right  in  Eq.  (70). 

According  to  the  observations  made  in  conjunction  with 
Eq.  (32),  the  occurrence  of  such  a  singularity  must  be  expected. 
One  deals  with  the  evaluation  of  a  certain  double  integral,  and 
such  a  singularity  will  arise,  unless  one  carries  out  the 
integrations  in  a  certain  sequence.  In  the  present  approach  a 
different  sequence  has  been  applied.  To  give  meaning  to  the 
double  integral  in  such  a  case,  one  first  excludes  the  point 
(x»y).  by  a  line  5=x  (1-e)  or  £=1-e>0.  In  the  present  context 

y0=°. 

The  factor  1 /p  in  the  integrand  arises  from  l_2,  Eq .  ( BM) . 
Before  the  limits  of  integration  are  substituted  the  term  which 
causes  this  factor  reads 


,  9  .  ,2  2.1  /2 
p-2(i-pV  ia_;p> — 

q-p 


(Ci) 


The  variable  of ^integration  q  is  expressed  by  £  in  Eq .  (B2).  The 
lower  limit  in  E,  in  Eq.  ( B1 )  ,  transforms  into  q=1  .  This  appears 
here  as  the  upper  limit.  The  upper  limit  is  q  =  |p|  for  values 
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of  p,  for  which  that  the  integration  over  5  ends  at  one  of  the 
straight  lines  AB  or  CB  in  Figure  14.  Along  the  straight  line 
BC,  the  upper  limit  is  given  by  by  5«1-e.  Then  one  obtains 

q  -  1-(1-e)-(1-p2) 

2  2 
q  =  e  +  p  -  ep  . 

A 

By  the  transition  ftom  5  to  q  the  upper  and  lower  limits  are 
interchanged . 

e  is  assumed  to  be  small,  and  obviously  p»0(e)  along 
BC.  The  last  term  on  the  right  can,  therefore,  be  neglected. 
The  values  of  p  for  which  this  limit  must  be  applied  are  found 
from  the  requirement,  that  at  the  transition  value  the  limits 
from  the  adjacent  regions  are  the  same 

1 P  |  -  e  +P2 

For  e  small  p2<<e. 

The  transition  points,  therefore,  lie  at 

I  P I  *  e- 

Substituting  these  limits  into  the  expression  (Cl)  one  obtains 
p“2(1-p2)"3/2  for  p  >  e 

?  i  /  2  2.1/2 

and  P  (1-p  )  [-- — 27172  - ' - - ]  f0r  P<e 

( 1-p  )  e-p 

In  the  bracket  we  have  used  the  fact  that  p  and  e  are  small. 
Notice,  however,  that  p/e  is  not  necessarily  small.  One  has 
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(e2-p2)'/2 


e-p 


_  1 _ r  2  ,  2  2.1/2i 

n  2. 172^  27[e_P  “(e  “P  J  J 

(1-p  )  (e-p  ) 

2  2  2  2 


1C*  C*  C*  Cm  M  - 

re  -2ep  +p  -e  +  p  2nf  h 

[  X  ■ 5  3.1  yo  P  01  oJ 


e  ‘  _2. /_2  2 »"l  /2 

e-p  +(e  -p  ) 


The  critical  integral  then  appears  in  the  forms 


-e  +e 

^  [ (a1-p)(a2+p)]  1/2fm,n(p)dp  +  |  p0(— dp 


1 

*  {  ^[(arp)(I2.p)]-1/2fn'n(p)<lp 


The  first  and  third  term  combined  are  the  integral 
interpreted  as  its  Cauchy  principal  value.  The  integrand  of  the 
second  term  can  be  split  into  its  symmetric  and  antisymmetric 
parts.  The  lowest  order  term  of  the  symmetric  part  is 


const 


0(e) 


This  establishes  the  (expected)  result. 
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APPENDIX  D 


EVALUATION  OF  AN  INTEGRAL 


The  expression  to  be  evaluated  is 


i 

P  |  [(a^pMag+p)]"1^'^?  ,  a.|  >0 ,  a2>0 

-a„ 


where  P  expresses  that  the  Cauchy  principal  value  is  to  be  taken, 


By  setting 


a1-a2  ~  a1+a2 

2  +  p  2 


one  obtains 


ai +a5  - 

i  ■  (-V1) 1  i 


where 


+  1 


-1 


dP. 


(1-p2)1/2(fS-a) 


and 


a  =  - 


ara2 

al+a2 


Be 


cause  a.|>0f  a2>0  one  has 


a  I  <  1 . 


Setting 


u2-1 
u2  +  1 


dp 


4udu 
(u2  +  1  )2 
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One  obtains 


/ 1  2  v  Hu41 

(1-p  )  -  —n - 

u  +1 


J  *T 
*  11 


2du 


+  CO 


u  (1-a)-(1+a) 


1  f  2du  1  f  du 
J  u2-b2  '  1-a  )  u2-t>2 


(D1  ) 


2  1  +a 

where  b  -  >  0,  since  | a  | <  1 


I  is  rewritten 


I  - 


1  1  f 

Ci-aT  2b  J  '■u-b 


u+F^du 


Eq.  ( D1 )  shows  that  the  path  of  integration  can  be  closed  at 
infinity,  in  the  upper  half  of  the  complex  u-plane  Figure  15. 
The  singular  points  u  »  _+  b  are  excluded  by  small  circles.  The 
integral  over  this  closed  path  gives  zero.  The  residue  at  the 
two  singular  points  are  _+1 .  The  contributions  of  the.  small 
circles  must  be  subtracted,  if  one  forms  the  principal  value  of 
the  integral.  The  sura  of  these  contributions  is  zero,  because 
the  residues  have  opposite  sign.  Therefore 


l1 


P  j  [ (a1 -p) (a2+p] j"1 /2p_1dp 


n 


r&rv  inKmirevyuix  w  >\  r 


APPENDIX  E 

LEIBNITZ'  RULE  FOR  CERTAIN  INTEGRALS  WITH  SINGULARITIES  IN 

THE  INTEGRAND 


Consider  an  integral 
x2(z) 

I(z)  -  |  f(x,z)dx 
x^z) 

where  the  function  f(x,z)  and  the  limits  x.|(z)  and  x2(z)  are 
differentiable  with  respect  to  z.  Then  one  obtains  according  to 
Leibnitz’  rule 

dl/dz  -  f (x2(z) ,z) (dx2/dz)  -  f (x- (z) ,z)(dx^/dz) 
x2(z) 

♦  |  Of  (x  ,z)/3z)dx  (El) 

X1  (z) 

To  derive  this  equation  one  introduces 
x  -  g(u,z) 

where  the  function  g(u,z)  is  chosen  in  such  a  manner  that  for  the 
values  of  z  under  consideration 

x1  (z)  -  g(u1  ,z) 

x2(z)  =  g(u2,z) 

where  u1  and  u2  are  independent  of  z. 

We  denote  by  f(1^  and  f(2^  or  g(1)  and  g(2)  the  derivative 
of  f  and  g  with  respect  to  either  the  first,  or  second  argument. 
For  higher  derivatives  a  corresponding  notation  is  used,  for 
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v  nm  n 


frsrn-i  n  m  A  *  A  M  'TfVtL  £4  ‘4! 


instance  for  the  second  derivative  of  f  with  respect  to 

the  first  argument. 

Then 


I(z)  -  j  f  (g(u.z)g^  (u,z)du 
U1 

and 

)  (g(u  ,z)  ,z)g  ^  (u  ,z)g  ^ 1  *  (u  ,z  )+f  2(g(u  ,z)  ,z  )g^ 1  ^  (uz) 

+  f  (g(u  ,z)g^ 1  ’  ^  (u  , z  )  }du 

One  obtains  the  result  Eq.  (El)  by  observing  that  in  a  more 
conventional  rotation 

f(1)g(1)  *=  (3f/3x)(3g/3u)=>(df/du)| 

z  -  const 

Carrying  out  an  integration  by  parts  in  the  first  term  of  the 
integrand  one  obtains 

u2  u  2 

dl/dz  -  f(g(u,z),z)  g(2)(u,z)|  ♦  j{-f(g(u(z),z)g(1,2)(u,z)  + 

U1  u  1 

+  f  ^  2^  (g  (u  ,z  )  ,z  )g^  (uz)  +  f (g (u ,z )g^ 1 ’ 2^ (u  ,z)  }du 


dl/dz 


U2 

[lf(1 


U, 


The  first  and  the  third  term  of  the  integrand  cancel. 


We  return  to  the  original  variable  x 
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f (x1  ,z ) 


f  (g(ufz)  ,z)  | 


u=u. 


g(2)(u,z) 


dx^/dz 


u-u . 


and  analogously  for  U2*  One  then  obtains  Eq.  (El). 

This  approach  is  useful,  if  the  function  f  has  a  singularity 
at  one  of  the  limits,  and  if  one  has  to  form  higher  derivatives. 
Consider  a  simple  example: 


z 

I(z)  -  |  (z-x)1/2dx 
o 


( E2) 


The  result  is  obvious 


z 


I(z)  -  -(2/3)(z-x)3/2  |  -  (2/3)z3/2 

(E3) 

x»o 

-  z1/2  ;  d2I/dz2  -  (1/2)z”1/2 

(E4) ,(E5) 

Let  us  forego  the  direct  integration,  and  carry  out  the 
differentiation  in  Eq.  (E2)  with  respect  to  z  according  to 
Leibnitz'  rule.  The  integrand  in  Eq.  (E2)  vanishes  at  the  upper 
limit  z. 

Accordingly, 


dl 

dz 


(1/2) 


(z-x)  1 ^2dx 


'(z-x) 


1  /  2 


1/2 


Although  the  integrand  in  Eq .  (E2)  has  a  square  root  singularity 
at  the  upper  limit,  one  obtains  the  correct  result,  Eq.  (E4). 

The  procedure  fails  for  the  second  derivative 
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z 


z 

(1/2)  (z-x)_1/2  j  -(1/4)  j  (z-x)~3/2dx 
o  o 


At  the  upper  limit  the  first  term  becomes  infinite.  It  is  true 
the  integral  becomes  infinite,  too,  and  the  two  " inf  ini  ties”  will 
cancel,  but  this  observation  still  does  not  give  a  result. 

Proceeding  in  the  manner  shown  above  we  introduce 


(x-z)  -  u 


Then 


z 

I(z)  ■  J  u1/2du 
o 

Because  the  integrand  has  no  singularity  at  the  upper  limit, 
Leibnitz*  rule  can  be  applied 


dl/dz  -  u1/2  j  -  z1/2 
o 

The  second  derivative  is  then  trivial. 

In  the  context  of  this  report,  the  problem  appears  in  a  more 
complicated  form. 

I(z2)  -  |  (x-a(z2)  J1 /2f  (x  ,z2)dx 
a  ( z  2 ) 

where  the  derivatives  of  f  and  a  with  respect  to  the  two 
variables  exist.  We  set 

2 

s  =  z 
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Then 


b 

I(s)  -  j  (x-a(s) J1 /2f (x  ,s)dx 
a(s) 

The  integrand  vanishes  at  the  lower  limit.  Application  of 
Leibnitz'  rule  gives  the  correct  result. 


b 

dl/ds  -  |  [-(1/2) (x-a(s))'1/2a'(s)f (x  ,s)  +  (x- (a(s)  ^  /2f(2) (x,s)]dx 
a  (s ) 

(E5) 


The  procedure  fails  for  the  second  derivative.  Therefore,  one 
sets 


x  -  a(s)  -  u 


Then 


1(3) 


-b-a(s) 

f 

=  J 

u =0 


u1 /2f (a(s)+u,s)du 


The  integrand  is  regular  at  the  upper  (s-dependent )  limit. 
Therefore  Leibnitz'  rule  can  be  applied 


1  /? 

(dl/ds)  =*  -u  f (a(s)  +u , s )a ' (s ) 

w=b-a(s) 


(E6) 


b-a(s) 

+  j  u1 /2[ f ( 1 ) (a(s)+u ,s)a' (s)  +  f ^ (a(s)+u,s]du 
u*=0 


To  obtain  Eq .  (E5)  an  integration  by  parts  must  be  carried  out. 
Notice  that 


(a(s)+u ,s)du 


f  (a(s)+u,s) 
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An  integration  by  parts  of  the  first  term  in  the  integrand  gives 
outside  of  the  integral  a  term  which  cancels  the  first  term  on 
the  right  in  Eq.  (E6).  Therefore 

b-a(s) 

*  |  [(1/2)u  ^2f  (a  (s  )+u  ,  s  )a ' (s )  +  ^2f  ^ (a(s)+a  ,s)  ]du 

u-0 

Eq .  (E5)  is  obtained  after  one  reintroduces  x.  The  second 
2  2 

derivative  (d  I/ds  )  from  Eq.  (E6)  can  be  found  without 
encountering  an  infinity. 

In  this  report,  one  needs  derivatives  with  respect  to  z  but 
only  for  z-0.  One  has 

(dl/dz)  -  2z(dl/du) 

(d2I/dz2)  -  2(dl/du)  +  4z2(d2I/du2) 


Hence 


dl/dz j  -  0 
z=0 


(E7) 


d2I/dz2  j  -  2(dl/du) 
z=0 


b 

-  {  -(x-a(z2)"1/2(da/dz2)f(x,(z)2)+2(x-a(t2)1/2f(2)(x,z2)dxl 

a (z2)  z" 


b 

a(z2) 


d 

d  (z2) 


[(x-(a(z2)1//2f(x,z2)]dx 

z  =0 


(E8; 


The  results  Eqs.  (E7)  and  (E8)  could  have  been  anticipated 
on  intuitive  grounds. 


80 


REFERENCE 


W.  Grobner,  M.  Hofreiter,  N.  Hofreiter,  J.  Laub,  and  E. 

II 

Peschl.  "Integraltaf el ,  Unbestimmle  Integrale,  Zentrale  fur 
wissenschaf tliches  Berlchtswesen  der  Luf tfahrtforschung 
Berlin-Adlershof  1944  (has  also  appeared  in  book  form). 


81 


Mach  Wave 


Figure  1 


(x,y) 


Region  of  Integration  with  Subsonic  and  Supersonic 
Leading  Edges. 
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Figure  3.  Subdivision  of  the  Region  of  Integration. 
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Figure  5.  Evaluation  of  the  Double  Integral  in  Eq,  (32)  in  the 
Vicinity  of  a  Control  Point. 
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Figure  6.  Region  of  Integration  for  z=0  with  Line  n  =  const 
Tangent  to  one  of  the  Leading  Edges. 
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Region  of  Integration  in  the  £ti  and  £  n  Coordinates 


Figure  10.  Subdivision  of  the  Region  II  if  the  Image  of  the 
Leading  Edge  has  a  Point  where  it  is  Tangent  to  a 
Line  n  =  const. 
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z^o  n=l 


♦ 

Figure  11.  Subdivision  of  Region  II  in  the  Cn  Plane  with  Further 
Subdivisions  in  the  £n  Plane. 
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Figure  12.  Regions  A  and  B  in  a  Problem  with  a  Subsonic  and  a 
Supersonic  Leading  Edge. 
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u 
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-b 
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+b 
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Figure  15.  Evaluation  of  a  Certain  Integral  in  the  Complex 
u^Plane. 
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